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SUMMARY 


Measurements  were  made  of  the  unsteady  normal  force  and  pitching 
moment  on  an  NACA  0012  airfoil  model  oscillated  both  sinusoidally  and 
nonsinusoidally  over  a  range  of  incidence  angles,  including  a  substantial 
penetration  into  stall.  The  sinusoidal  normal  force  and  pitching  moment 
data  were  reduced  and  tabulated  as  functions  of  the  angle  of  attack,  the 
angular  velocity  parameter,  and  the  angular  acceleration  parameter.  This 
generalized  form  of  the  data  was  used  to  reconstruct  the  measured  sinu¬ 
soidal  aerodynamic  response  of  the  model  airfoil  with  excellent  results. 
Additional  correlations  were  made  using  nonsinusoidal  pitch  schedules 
which  included  periodic  ramp  changes  in  angle  of  attack  and  a  flexured 
angular  blade  response  to  a  one-per-rev  sinusoidal  incidence  angle  change 
typical  of  that  for  a  helicopter  blade.  The  agreement  between  predicted 
and  measured  normal  force  and  moment  loops  was  very  good  for  the  ramp 
motion.  Integration  of  the  nonsinusoidal  moment  loops  to  calculate  aero¬ 
dynamic  damping  in  pitch  has  produced  good  agreement  between  predicted  and 
measured  results.  Good  agreement  was  also  obtained  between  predicted  and 
measured  flexible  blade  response. 

Scaling  procedures  describing  aspect  ratio  and  compressibility  effects 
were  derived  and  used  to  calculate  full-scale  unsteady  normal  force  and 
moment  coefficients  from  the  unsteady  model  data.  Correlation  studies 
were  performed  to  determine  whether  the  incorporation  of  the  unsteady 
aerodynamic  data  into  an  existing  rotor  aeroelastic  analysis  would  provide 
an  adequate  method  for  predicting  stall -flutter -induced  control  loads. 

The  analysis  was  carried  out  for  the  NH-3A  (S-61F)  and  CH-53A  helicopter 
rotors  and  for  a  wind  tunnel  model  rotor .  Use  of  variable  rotor  inflow 
model  and  unsteady  aerodynamics  provided  good  correlation  with  torsional 
oscillations  observed  in  test. 
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This  work  was  performed  for  the  Eustis  Directorate,  U.  S.  Army  Air 
Mobility  Research  and  Development  Laboratory,  Fort  Eustis,  Virginia, 
under  Contract  DAAJ02-71-C-0003 ,  DA  Task  1F162204A13904 .  Technical 
monitor  for  the  Army  was  Mr.  William  Nettles.  The  experimental  work  and 
the  unsteady  airfoil  data  analysis  were  performed  at  the  United  Aircraft 
Research  Laboratories,  East  Hartford,  Connecticut,  by  Messrs.  F.  0.  Carta 
and  G.  L.  Commerford.  The  helicopter  rotor  control  load  correlations 
were  performed  at  the  Sikorsky  Aircraft  Division  of  United  Aircraft 
Corporation,  Stratford,  Connecticut,  by  Dr.  R.  G.  Carlson  and 
Mr.  R.  H.  Blackwell. 
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Classical  rotor  theories  are  incapable  of  accurately  predicting 
overall  rotor  performance  characteristics  in  stall  (References  1  through 
7),  nor  can  they  predict  the  details  of  the  rotor  loading  required  for  the 
definition  of  critical  high-frequency  vibratory  loads.  Such  theories 
generally  assume  that  two-dimensional,  steady-state  airfoil  data  are 
applicable  and,  further,  that  the  inflow  velocities  induced  at  the  rotor 
blades  by  the  vortex  system  of  the  rotor  are  constant  with  respect  to  time 
and  space.  In  an  attempt  to  provide  a  more  suitable  analysis  for  applica¬ 
tion  to  the  rotor  stall  situation,  a  study  was  performed  at  the  United 
Aircraft  Corporation  (UAC),  under  USAAMRDL  sponsorship,  to  integrate 
available  two-dimensional,  unsteady  airfoil  data  into  an  existing  rotor 
aeroelastic -variable  inflow  analysis  (Reference  8).  The  data  available 
for  this  purpose  were  those  obtained  from  Reference  9  for  airfoils 
executing  pure  sinusoidal  pitch  motion.  Since  rotor  blades  do  not 
normally  execute  pure  sinusoidal  motion  as  they  penetrate  the  tall  regime, 
a  unique  feature  of  the  expanded  UAC  analysis  was  the  conversion  of  the 
available  sinusoidal  results  to  a  more  general  nonsinusoidal  form  in  which 
local  section  forces  and  moments  are  expressed  in  terms  of  instantaneous 
Mach  number,  angle  of  attack,  a,  angular  velocity,  A,  and  angular  acceler¬ 
ation  parameter,  B.  This  conversion  is  required  before  section  forces  and 
moments  associated  with  typical  rotor  blade  angle -of -attack  variations  can 
be  computed  in  a  rational  manner.  Figure  1  shows  a  typical  angle-of -attack 
variation  for  one  rotor  blade  station  along  with  the  associated  variation 
in  angular  velocity  parameter.  The  angular  velocity  parameter  has  been 
shown  (References  10  and  11,  for  example)  to  have  an  important  influence 
on  airfoil  stall  characteristics.  The  departure  from  pure  sinusoidal 
motion  is  evident,  particularly  in  the  critical  angular  velocity  parameter. 

Reference  7  contained  data  similar  to  those  of  Reference  9  but  at 
higher  subsonic  Mach  numbers.  Attempts  were  made  to  recast  these  com¬ 
pressible  data  in  a  more  general  form  dependent  on  a,  A,  and  B.  However, 
the  range  of  parameters  tested  in  Reference  7  was  too  limited  to  permit 
this  transformation.  For  this  reason,  the  UAC  rotor  analysis  developed  in 
Reference  8  used  the  incompressible  data  of  Reference  9  in  conjunction 
with  scaling  laws  to  account  approximately  for  compressibility.  Although 
the  accuracy  of  the  analysis  was  compromised  to  a  certain  degree,  its 
ability  to  predict  a  stall  flutter  type  of  response  was  successfully 
demonstrated.  However,  correlation  of  analytical  and  flight  test  data  on 
the  NH-3A  helicopter  indicated  that  the  predicted  stall  flutter  response 
occurred  at  rotor  lifts  higher  than  those  measured  in  flight.  Although 
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this  correlation  with  NH-3A  flight  data  was  not  good,  a  comparable 
correlation,  using  the  same  a,  A,  B  technique,  with  model  scale  rotor  data 
yielded  excellent  agreement  (see  Figure  12  of  Reference  12). 

To  achieve  good  correlation  between  the  predicted  response  from  the 
analysis  of  Reference  8  and  the  measured  response  on  a  full-scale  helicop¬ 
ter  rotor,  it  is  necessary  that  (l)  reliable  detailed  flight  test  data  be 
available,  (2)  all  facets  of  the  aeroelastic  analysis  (e.g.,  prediction  of 
rotor  inflow,  applicability  of  the  sinusoidal  unsteady  data,  scaling  of 
data  to  account  for  compressibility  effects,  representation  of  blade  flex¬ 
ibility,  etc.)  be  accurate,  and  (3)  factors  neglected  in  the  correlation 
(e.g.,  effects  of  blade  three-dimensional  flow,  wake  geometry  distortions, 
or  fuselage /wing  interference)  be,  in  fact,  negligible. 

It  was  decided  to  examine  the  effects  of  each  of  these  factors  singly 
if  possible.  A  logical  starting  point  for  this  examination  concerned  the 
fundamental  assumption  that  unsteady  sinusoidal  data,  tabulated  as  a 
function  of  a,  A,  and  B,  can  be  used  to  predict  the  aerodynamic  response 
of  an  airfoil  executing  a  nonsinusoidal  motion.  To  test  this  assumption, 
a  group  of  proof -of -theory  experiments  was  devised  and  carried  out  as 
reported  herein.  Concurrently,  an  analysis  was  carried  out  in  an  attempt 
to  improve  the  correlation  of  predicted  and  measured  control  loads.  The 
new  unsteady  data  tables  were  incorporated  into  the  UAC  aeroelastic  blade 
response  program,  and  correlation  studies  were  carried  out  for  the  NH-3A 
and  the  CH-53A  rotor  systems  and  for  a  dynamically  scaled  model  rotor. 

The  effect  of  rotor  inflow  upon  torsional  response  prediction  was  assessed 
by  repeating  the  analysis  for  several  inflow  distributions.  Slight  modi¬ 
fications  were  made  to  the  unsteady  data  tabulations  of  Reference  8  to 
smooth  out  regions  of  irregular  data.  This  modified  original  unsteady  data 
and  the  new  model  data  were  then  compared  by  running  some  of  the  correla¬ 
tion  cases  with  both  tabulations. 


MODEL  AIRFOIL  AND  TEST  FACILITY 

As  stated  before,  one  of  the  primary  objectives  of  this  experimental 
program  was  to  test  the  validity  of  the  hypothesis  that  sinusoidal  data 
could  be  generalized  by  the  UAC  a,  A,  B,  method  (References  8  and  9)  to 
predict  the  aerodynamic  response  of  an  airfoil  to  a  nonsinusoidal  motion. 
This  was  designed  to  be  a  completely  self-contained  program  in  which  a 
single  airfoil  would  be  tested  in  the  same  facility  with  both  sinusoidal 
and  nonsinusoidal  motions  to  eliminate  as  many  extraneous  variables  as 
possible  in  the  course  of  this  proof -of -theory  test. 

In  view  of  the  multitude  of  tasks  required  of  the  model  airfoil,  there 
were  a  number  of  design  constraints  imposed  on  its  internal  and  external 
configuration.  First,  it  was  required  that  the  test  program  be  performed 
in  a  convenient  test  facility  in  which  a  large  number  of  varied  tests 
could  be  run  over  a  long  period  of  time  without  incurring  occupancy  con¬ 
flicts.  The  facility  chosen  to  satisfy  these  requirements  was  the  UAC 
two-dimensional,  high-speed  cascade  wind  tunnel  (shown  in  Figure  2)  modi¬ 
fied  substantially  for  the  performance  of  this  isolated  airfoil  test. 

This  tunnel  received  air  from  a  large  pressurized  plenum  chamber.  After 
passing  through  a  series  of  fine  mesh  screens  to  reduce  turbulence,  the 
air  entered  the  test  section  through  a  rectangular  bellmouth.  The  test 
section  was  4  inches  wide,  and  the  adjustable  sidewall  quadrants  were 
rotated  to  provide  a  test  section  height  of  14  inches  between  the  horizon¬ 
tal  floor  and  ceiling  pieces.  Sidewall  extension  pieces  were  added  to  the 
tunnel,  and  the  air  exhausted  to  atmosphere  four  chordlengths  downstream 
of  the  model  airfoil  trailing -edge  location. 

The  specific  tasks  required  of  the  airfoil  model  each  placed 
additional  constraints  on  the  model  design.  The  maximum  4 -inch  width  of 
the  tunnel  and  the  need  for  sidewall  boundary  layer  fences  (see  below) 
restricted  the  total  model  span  to  3-73  inches.  The  chord  of  the  model 
was  chosen  to  oe  5  inches  to  permit  the  insertion  of  10  miniature  pressure 
transducers.  The  need  for  pressure  transducers  located  along  the  entire 
chord  produced  a  high  strength  requirement  along  the  entire  chord,  but  the 
need  for  a  flexured  oscillating  airfoil  scaled  to  simulate  a  full-scale 
rotor  blade  section  produced  a  low  inertia  requirement.  Accordingly,  the 
model  was  constructed  as  described  below  to  fulfill  as  many  of  these 
requirements  as  possible. 


The  basic  building  blocks  of  this  model  were  the  steel  spar  and  the 
magnesium  ribs,  shown  partially  assembled  in  Figure  3*  The  large  magnesium 
rib  contained  all  of  the  pressure  instrumentation.  In  the  upper  half  of 
Figure  3,  this  portion  of  the  rib  is  seen  to  contain  10  transverse  holes 
into  which  the  10  miniature  pressure  transducers  were  inserted.  In 
addition,  two  alignment  and  drive  pins  are  seen  which  were  matched  to  holes 
on  the  narrow  rib;  this  constituted  the  driving  half  of  the  model.  The 
lower  half  of  Figure  3  shows  the  back  side  of  the  large  rib.  which  also 
contained  20  static  pressure  taps,  10  for  each  surface.  Also  shown  here 
is  the  steel  spar,  hollowed  out  to  compensate  for  the  model  weight  aft  of 
the  25  percent  chord  pivot  axis.  Instrumentation  leads  and  hypo  tubing 
were  passed  through  the  hollow  center  of  the  spar,  and  the  round  shaft  ends 
were  later  inserted  into  support  bearings  in  both  sidewalls . 

Balsa  filler  pieces  were  then  placed  fore  and  aft  of  the  shaft  as 
shown  in  Figure  4.  Balsa  wood  was  specifically  chosen  to  provide  a  light¬ 
weight  rigid  core.  Next,  the  two  model  halves  were  assembled  and  four 
layers  of  carbon-epoxy  skin  were  wrapped  around  the  airfoil :  two  with 
fibers  in  the  chordwise  direction,  and  two  with  fibers  in  the  spanwise 
direction.  The  static  pressure  tubing  which  can  be  seen  protruding  from 
the  surface  of  the  rib  in  Figure  4  was  partially  buried  in  the  skin  and 
then  machined  flush  with  the  surface.  After  the  surface  was  smoothed,  the 
skin  was  sliced  chordwise  at  the  joint  between  the  two  ribs,  the  mating 
faces  were  machined  flush,  and  the  model  was  painted.  Static  holes  were 
then  drilled  from  both  airfoil  surfaces  into  the  transducer  cavities. 

This  is  shown  in  Figure  5. 

The  miniature  differential  pressure  transducers  were  supplied  by 
Sensotec  (model  no.  SA-SD,  Type  8j-6h)  and  had  a  working  range  of  *15  psi. 
Their  sensitivity  at  unit  gain  varied  from  1  mv/psi  to  2.3  mv/psi,  and 
the  actual  values  as  calibrated  are  given  in  Table  I.  Although  the  trans¬ 
ducers  showed  some  deviation  with  temperature  charge,  these  deviations 
were  small,  and  the  characteristic  was  linear  at  constant  temperature. 

This  is  shown  in  Figure  6  for  a  typical  transducer  at  a  gain  factor  of  20. 
It  will  be  shown  later  that  the  mode  of  operation  chosen  avoided  all 
possible  thermal  problems. 

Ore  of  the  transducers  is  pictured  in  the  upper  half  of  Figure  7. 

The  transducer  consisted  of  a  circular  metal  diaphragm  l/8  inch  in  diam¬ 
eter,  with  bonded  semiconductor  strain  gages.  The  diaphragm  was  mounted 
at  one  end  of  a  cylindrical  tube,  which  was  partially  open  at  the  far  end 
and  capped  by  an  overhanging  end  piece.  An  0-ring  on  the  cylinder  served 
to  isolate  the  two  ends  of  the  transducer  when  it  was  inserted  into  its 
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c  Wty  in  the  model.  A  typical  installation  is  shown  schematically  in  the 
lcwer  half  of  Figure  7.  Pictured  here  is  a  sketch  of  a  spanwise  section 
through  the  large  rib,  looking  aft.  The  figure  is  largely  self-explanatory. 
The  epoxy  plugs  at  the  two  surfaces  were  poured  in  place  to  prevent 
pressure  leakage  through  the  layers  of  skin.  All  but  the  last  three  trans¬ 
ducers  were  configured  in  this  manner.  Because  of  space  limitations,  the 
three  rearmost  transducers  were  placed  forward  of  their  measuring  stations 
and  were  connected  to  them  by  means  of  short  surface  tunnels  buried  under 
the  skin.  These  are  seen  in  Figures  3  and  4  as  short  chordwise  streaks  in 
the  large  rib.  Both  the  transducer  and  static  tap  orifices  were  distribu¬ 
ted  along  the  airfoil  chord  according  to  the  normalized  roots  of  a 
Legendre  polynomial  to  take  advantage  of  the  inherent  accuracy  of  the 
Gaussian  integration  method  (References  13  and  14).  This  is  discussed  in 
a  later  section  and  in  Appendix  I. 

Sidewall  boundarv  layer  end  plates  were  used  in  an  attempt  to  reduce 
sidewall  interaction  effects  on  such  a  short  span  model.  The  end  plates 
were  9.25  inches  long  and  3  inches  high,  with  a  circular  leading  edge,  as 
shown  in  Figure  8.  They  were  located  0.188  inches  from  the  sidewalls, 
with  their  leading  edges  2  inches  forward  of  the  airfoil  leading  edge. 
Installation  of  the  end  plates  produced  an  effective  aspect  ratio,  based 
on  reduction  of  lift  curve  slope,  of  approximately  9.5  (see  page  84  of 
Reference  15). 

The  airfoil  model  was  mounted  in  the  wind  tunnel  between  thrust 
bearings  which  were  loaded  to  keep  the  two  airfoil  halves  together.  The 
airfoil  was  oscillated  either  sinusoidally  or  nonsinusoidally  by  the  drive 
system  shown  schematically  in  Figure  9.  In  this  system,  a  variable -speed 
6.5-hp  electric  motor  drove  an  eccentric  cam,  located  in  a  transmission 
box,  in  a  circular  motion  by  means  of  a  drive  belt  and  pulleys.  A  double 
follower  within  the  transmission  converted  the  circular  cam  motion  to  a 
periodic  shaft  motion,  which  was  transmitted  to  the  blade  through  a  shaft 
coupler.  Tor  mo.,t  of  the  tests,  this  coupler  was  rigid  and  the  airfoil 
motion  was  completely  prescribed  by  the  angular  schedule  of  the  cam- 
follower  system.  For  a  small  portion  of  the  test  program,  the  coupler 
was  flexible.  In  these  cases,  the  model  was  driven  at  low  frequency 
(e.g.,  approximately  10  cps)  to  simulate  a  one -per -rev  cyclic  change  in 
angle  of  attack  of  a  typical  rotor  blade,  and  the  coupler  flexibility  was 
chosen  to  simulate  the  torsional  response  of  a  typical  rotor  blade  at 
frequencies  between  three -per -rev  and  ten-;»er-rev.  Most  of  the  tests  were 
conducted  with  a  sinusoidal  cam  which  imparted  an  angular  amplitude  of 
*8  degree  to  the  blade.  For  a  portion  of  these  tests,  the  pitching 
oscillation  either  approximated  the  time  history  of  a  helicopter  rotor 
blade  or  executed  a  ramp-like  motion,  both  at  *8  degrees.  A  very  small 
number  of  sinusoidal  tests  were  conducted  at  an  amplitude  of  *3  degrees. 
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Further  descriptions  of  each  test  condition  are  given  later  in  the  text, 
and  an  overall  summary  of  the  test  points  actually  run  is  given  in 
Table  II. 


Angular  deflections  were  measured  using  a  cylindrical  electrical 
capacitor  on  each  side  of  the  airfoil.  One  part  was  mounted  on  the  air¬ 
foil  shaft  and  rotated  with  the  same  motion  as  the  shaft.  The  other  part 
was  stationary  and  its  cylindrical  surface  was  made  concentric  with  the 
shaft  center  line.  This  is  shown  schematically  in  Figure  10.  Both 
cylindrical  surfaces  were  coated  with  conducting  paint  and  were  excited  by 
means  of  a  capacitive  signal  conditioner.  An  initial  setting  was  chosen 
with  the  faces  staggered  by  some  amount,  as  shown  in  the  sketch  of  Figure 
10.  As  the  shaft  rotated  in  either  direction,  there  would  be  more  or  less 
capacitive  area  in  opposition,  which  was  manifested  as  a  variable  output 
voltage.  This  voltage  was  found  to  be  a  linear  function  of  the  angular 
displacement  from  rest  with  sufficient  sensitivity  to  ensure  a  reading 
error  of  less  than  0.1  degree. 

One  of  these  capacitive  transducers  was  mounted  at  the  free  end  of 
the  blade  and  always  measured  blade  motion.  The  other  transducer  was 
mounted  between  the  transmission  and  the  shaft  coupler  and  always  measured 
the  motion  input  to  the  coupler.  For  a  rigid  coupling,  the  two  trans¬ 
ducers  yielded  the  3ame  angular  motions,  and  only  one  was  recorded.  For 
the  flexible  coupling  tests,  the  coupler  transducer  still  measured  the 
input  motion,  but  the  blade  transducer  measured  the  effects  of  blade 
inertia,  blade  stiffness,  and  aerodynamic  forces  and  moments.  Hence,  the 
blade  transducer  produced  a  response  that  was  measurably  different  from 
that  of  the  coupler  transducer  during  a  flexible  test,  and  both  responses 
were  recorded. 


The  instantaneous  angle  of  attack  a  at  time  t  during  a  sinusoidal 
oscillation  of  frequency  <w  rad/sec  can  be  expressed  as 


where  Oty  and  a  are  the  mean  angle  of  attack  and  the  angular  amplitude  of 
oscillation,  respectively.  If  simultaneous  time -varying  measurements  of 
a  and  Ap  from  the  chordwise  pressure  transducers  described  in  the 
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preceding  section  are  made,  the  time-varying  lift  and  moment!  versus  angle 
of  attack  can  be  obtained.  The  normal  force  coefficient  (positive  up)  at 
any  instant  is  given  by 


'  \  •  ' 

and  the  corresponding  moment  coefficient  about  art  axis  at  a  (positive  nose 

up)  is  given  by  1 


-°HPU-PL 


(3) 


where  pu  and  p^  are  the  static  pressures  on  the  upper  and  lower  surfaces, 
respectively  (Ap  =  pu  -  Pl)»  c  is  the  blade  chord,  and  q  is  the  free- 
stream  dynamic  pressure .  1 

i 

» 

The  initial  portion  of  this  program  was  concerned  with  the  formulation 
of  a  tabulation  of  unsteady  normal  force  and  moment  for  an  airfoil 
oscillating  sinusoidally  in  pitch.  The  normal  force  and  moment  were  to  be 
tabulated  as  functions  of  angular  position,  a,  angular  velocity  parameter1, 


i 


and  angular  acceleration  parameter, 

■  i  .  i 


as  the  independent  variables.  Using  Equations  (l)  through  (3),  it  was 
possible  to  obtain  unsteady  normal  force  and  moment  data  in  terms  of  the 
experimental  variables  a,  a a ,  and  u>.  As  shown  in  Reference  &,  the 
required  tabulation  can  be  constructed  from  these  data. 
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The  oscillatory  portions  of  the  unsteady  pressures  were  measured 
separately  from  steady  or  mean  pressures.  The  steady  pressures  were 
obtained  directly  from  the  20  static  surface  taps.  The  oscillatory 
pressures  were  obtained  from  the  10  differential  transducers  by  (a)  first 
setting  the  tunnel  speed  at  the  desired  value  with  the  model  airfoil  set 
at  its  mean  angle  of  attack,  (b)  electronically  nulling  the  voltage  output 
of  the  transducers,  and  (c)  oscillating  the  airfoil  and  recording  the 
oscillatory  pressures  only  (see  below  for  further  details).  This  procedure 
was  followed  for  two  reasons.  First,  as  shown  in  Figure  6,  the  trans¬ 
ducers  were  found  to  be  temperature  sensitive.  Although  this  sensitivity 
was  small,  it  could  have  led  to  measurement  errors  because  the  air  temper¬ 
ature  with  flow  was  different  from  the  ambient  temperature  of  the  room  air. 
By  nulling  the  signal  under  full  flow  conditions  with  the  tunnel  at  its 
equilibrium  temperature,  the  thermal  drift  problem  was  eliminated.  Second, 
by  nulling  out  the  steady  portion  of  the  signal,  it  was  possible  to  use 
the  entire  voltage  range  of  the  transducers  to  measure  the  oscillatory 
portion  of  the  unsteady  pressure,  which  increased  the  accuracy  and  sensi¬ 
tivity  of  these  measurements. 

The  steady-state  tests  were  run  with  the  airfoil  restrained  at  the 
transmission.  The  steady  pressure  difference  distribution  was  obtained 
from  the  static  pressure  taps.  Although  the  differential  pressure  trans¬ 
ducers  were  not  used  for  measuring  the  steady  portion  (mean  value)  of  the 
pressure  during  the  oscillatory  tests,  it  is  noteworthy  that  they  were 
indeed  capable  of  measuring  the  steady  portion  accurately  if  the  trans¬ 
ducers  were  initially  balanced  at  the  operating  temperature  of  the  tunnel. 
This  can  be  shown  by  comparing  a  typical  steady-state  pressure  difference 
distribution  obtained  from  the  pressure  transducers  (with  proper  thermal 
compensation  and  calibration)  to  that  obtained  directly  from  the  static 
pressure  taps.  These  data  are  shown  in  Figure  11. 

The  steady-state  normal  force  and  moment  at  the  angles  of  attack 
listed  in  Table  II  were  obtained  by  integrating  the  static  pressure 
distributions  over  the  airfoil.  The  details  and  results  of  this  integra¬ 
tion  will  be  covered  in  a  later  section.  The  pressures  were  recorded  both 
from  a  manometer  board  and  from  a  pressure  scanner  system.  The  pressure 
scanner  (scanivalve)  mode  of  operation  allowed  the  sequential  monitoring 
of  the  static  pressures.  Each  pressure  input  was  scanned  in  sequence  and 
converted  to  an  analog  voltage  with  a  single  pressure  transducer  external 
i  to  the  model.  The  voltage  output  was  displayed  on  a  digital  voltmeter. 


The  remainder  and  bulk  of  the  test  program  was  concerned  with 
obtaining  unsteady  normal  force  and  moment  data  on  the  airfoil  while  it 
was  undergoing  either  sinusoidal  or  nonsinusoidal  pitch  oscillations. 

These  different  motions  were  produced  by  changing  cams  in  the  transmission 
or  by  replacing  the  rigid  coupler  with  one  of  three  separate  Bendix 
flexures  (series  $000).  The  test  procedure  used  with  each  of  these 
changes  was  essentially  the  same  and  is  described  below. 

First,  the  airfoil  mean  angle  of  attack  was  set  to  some  desired 
value,  corresponding  to  the  mean  position  of  the  transmission  output  shaft. 
All  angles  were  measured  with  a  Model  B  Ailger  and  Watts  Vernier  Angle 
Gage  (Clinometer)  which  was  calibrated  to  within  one  minute  (1')  of  angular 
deviation.  A  special  block  was  machined  to  fit  closely  over  the  airfoil 
with  its  upper  surface  parallel  to  the  airfoil  center  line  (see  Figure  12). 
Angular  measurements  were  made  with  the  Clinometer  positioned  on  this 
block  as  shown. 

With  the  airfoil  positioned  as  described,  the  wind  tunnel  was  started. 
The  test  section  Mach  number  was  set  nominally  for  M  =  0.325  (actual 
velocity  at  operating  temperature  was  V  =  390  - 6  ft/sec).  The  tunnel  was 
run  long  enough  for  the  airfoil  to  reach  a  thermal  equilibrium.  At  this 
point,  the  tunnel  was  shut  down  quickly  and  all  the  transducers  were 
balanced  to  zero  output  with  temperatures  at  or  very  close  to  their 
operating  temperatures.  The  tunnel  was  restarted  and  the  Mach  number  was 
set.  This  procedure  eliminated  all  or  most  of  the  DC-bias  in  the  trans¬ 
ducer  output  due  to  temperature.  It  also  allowed  the  gages  to  be  balanced 
without  thermal  offset  so  that  the  output  signal  could  be  amplified  as 
much  as  possible  without  risk  of  saturating  the  recording  equipment  due  to 
thermal  drift.  The  transmission  locking  mechanism  was  released,  and  the 
motor  which  oscillated  the  airfoil  was  started.  The  required  frequencies 
were  set  by  monitoring  the  position  gage  output  on  a  Beckman  Universal 
Model  8370  EFUT  (Events  Per  Unit  Time)  Meter  and  adjusting  the  variable 
speed  control.  The  basic  system  used  for  recording  the  data  is  discussed 
in  Reference  16  and  is  referred  to  as  the  WISARD  (Wideband  System  for 
Acquiring  and  Recording  Data).  The  various  components  of  the  WISARD  and 
the  peripheral  equipment  used  throughout  this  program  are  shown  in  the 
block  diagram  of  Figure  13.  The  pressure  transducer  signals  were  ampli¬ 
fied  to  t  maximum  level  without  exceeding  l.Ov  output  (by  a  gain  of  50  near 
the  leading  edge  and  200  near  the  trailing  edge).  High-frequency  noise 
was  minimized  by  filtering  the  amplifier  output  above  10,000  cps.  Since 
this  frequency  was  two  orders  of  magnitude  higher  than  the  highest  oscilla¬ 
tory  test  frequency,  no  signal  attenuation  in  the  frequency  range  of 
interest  was  possible.  The  amplified,  filtered  pressure  transducer 
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signals  together  with  the  balanced  position  gage  signals  were  recorded 
directly  on  FM  tape  using  a  CEC  (model  VR-3400)  14  -channel  variable -speed 
tape  recorder.  Ten  channels  were  used  for  the  pressure  data,  and  two 
were  used  for  the  angle-of -attack  data  from  the  position  gages.  (The 
other  two  channels  were  used  for  identification  and  internal  mode  control 
of  the  WISARD.) 

Typical  raw  data  output  from  a  pressure  transducer  near  the  leading 
edge  (Ap)  and  the  position  gage  (a)  are  shown  in  Figure  14.  All  the 
panels  shown  in  this  figure  are  for  -8  degrees  sinusoidal  motion  at  a 
frequency  of  25  cps  and  a  Mach  number  of  0.325 .  The  left-hand  panels  are 
for  a  mean  angle  of  attack  of  6  degrees,  where  no  dynamic  stall  is 
observed.  The  right-hand  panels  are  for  a  mean  angle  of  attack  of  11 
degrees,  where  considerable  penetration  into  stall  can  be  observed.  Both 
lower  panels  are  oscilloscope  photographs  of  the  pressure  and  position 
time  histories  directly  from  the  tape  recorder  playback.  The  upper  panels 
are  Lissajous  figures  of  these  same  time  histories  in  which  the  position 
signal  was  used  as  the  horizontal  base.  The  characteristic  separation 
loop  is  evident  in  the  right-hand  panel. 

To  compute  the  unsteady  normal  force  and  moment,  the  raw  pressure 
data  were  multiplied  by  their  respective  sensitivities  and  divided  by  the 
gains  used  in  the  amplifiers.  The  mean  value  of  the  pressure  time 
histories  was  subtracted  from  the  actual  value  at  each  instant  of  time  so 
that  the  resulting  normal  force  and  moment  were  purely  the  unsteady  por¬ 
tions.  The  steady-state  normal  force  and  moment  for  a  given  angle  of 
attack  obtained  earlier  were  then  added  to  these  unsteady  data  to  obtain 
the  total  normal  force  and  moment . 

The  data  processing  procedure  required  that  analog  signals  from  the 
FM  tape  be  converted  to  digital  form.  This  conversion  was  also  done  on 
the  WISARD  and  is  shown  schematically  in  the  lower  portion  of  the  block 
diagram  of  Figure  13.  Here  the  operating  flexibility  of  the  data  system 
permitted  us  to  maximize  the  frequency  resolution  of  each  unsteady  signal, 
regardless  of  the  oscillatory  frequency.  By  recording  the  data  at  any  of 
a  number  of  selected  tape  speeds  and  playing  the  tape  back  at  a  slower 
speed,  it  was  possible  to  obtain  test  information  at  higher  frequencies 
than  would  otherwise  be  possible  if  recording  and  playback  speeds  were 
equal.  Once  the  recording  speed  was  selected,  the  amount  of  record  which 
would  fill  the  computer  storage  could  be  established  directly  from  the 
samping  rate  of  the  WISARD  digitizing  process.  A  detailed  description  of 
the  logic  used  in  selecting  tape  recorder  speeds  and  record  lengths  is 
given  in  Appendix  II.  The  results  of  this  study  are  shown  in  Table  III. 


These  settings  were  established  so  that  if  a  Fourier  analysis  of  the 
signals  was  required,  there  would  be  at  least  four  samples/cycle  in  the 
tenth  harmonic  of  the  oscillatory  frequency. 


TEST  RESULTS 


DATA  PROCESSING 


The  digitized  magnetic  tape  produced  by  the  data  system  was  processed 
on  a  UNIVAC  1108  digital  computer.  First,  the  millivolt  equivalents  of 
pressure  difference,  Ap,  and  angular  position,  a,  were  arranged  in  the 
computer  core  as  functions  of  time.  The  pressures,  suitably  multiplied  by 
conversion  factors,  were  integrated  to  yield  normal  force  and  moment 
coefficient  time  histories,  Cn(t)  and  Cm(t),  using  the  nondimensional 
versions  of  Equations  (2)  and  (3).  Because  the  pressure  orifices  were 
arranged  in  a  Gaussian  net  over  the  forward  90  percent  of  the  blade  chord, 
it  was  possible  to  take  advantage  of  the  Gaussian  integration  method,  as 
discussed  in  Appendix  I.  By  this  means,  a  greater  degree  of  accuracy  was 
attained  than  would  be  possible  using  standard  numerical  quadrature  formu¬ 
las.  Time  histories  representing  50  or  more  cycles  of  motion  were  con¬ 
verted  and  stored  for  each  test  point. 

These  data  contained  some  noncoherent  noise  in  addition  to  the 
desired  unsteady  information,  particularly  during  flow  separation.  A 
typical  time  history  for  separated  flow  is  shown  in  Figure  15.  To  minimize 
the  effects  of  this  superimposed  noise  and  to  enhance  the  coherent  or 
repeatable  elements  of  the  time  history,  a  Uo-cycle  signal  averaging  pro¬ 
cess  was  used.  This  reduced  the  noise  level  relative  to  the  desired  sig¬ 
nal  level  by  a  factor  of -^40  ~  6 .3  (see  pages  116-118  of  Reference  17). 

The  resulting  smoothed  cycle  for  each  component  was  then  normalized  with 
respect  to  its  period  to  occupy  100  units  of  nondimensional  time  to  provide 
a  uniform  data  set  for  further  processing.  A  sample  of  the  smoothed  time 
histories  is  shown  in  Figure  16  for  the  same  conditions  that  produced 
Figure  15. 

The  normal  force  and  moment  time  histories  for  each  test  point  were 
then  interpolated  at  integral  values  of  incidence  angle,  and  the  results 
were  stored  in  the  computer.  From  this  point,  the  process  described  in 
Reference  8  to  convert  the  normal  force  and  moment  into  tabular  functions 
of  angular  position,  a,  angular  velocity  parameter,  A,  and  angular 
acceleration  parameter,  B,  was  followed.  For  convenience,  the  schematic 
outline  of  the  process,  from  Figure  70  of  Reference  8,  is  reproduced  herein 
as  Figure  17.  Reference  8  describes  the  conversion  process  in  detail,  and 
this  description  will  not  be  repeated  herein.  The  generalized  unsteady 
normal  force  coefficient  r,  contained  in  Table  IV,  and  the  moment  coeffi¬ 
cient  is  in  Table  7.  Within  each  table  the  data  are  arranged  in  blocks 
for  each  integral  value  of  ge  metric  angle  of  ittack. 


STEADY-STATE  NORMAL  FORCE  AND  MOMENT 


The  oscillatory  tests  were  designed  to  provide  only  the  unsteady 
portions  of  the  normal  force  and  moment.  Therefore,  it  was  important 
throughout  this  program  that  an  accurate  measurement  of  the  steady  portion 
of  the  force  and  moment  be  available.  This  was  particularly  true  because 
of  the  small  span  of  the  model,  which  made  it  impossible  to  use  results  of 
previous  investigators.  Furthermore,  at  high  angles  of  attack,  the  steady 
portions  of  the  total  signals  were  considerably  greater  than  the  unsteady 
portions . 

To  ensure  the  accuracy  of  the  steady  force  and  moment,  these  data 
were  checked  in  as  many  ways  as  possible.  First,  the  pressure  measurements 
were  made  at  chordwise  grid  locations  corresponding  to  a  modified  Gaussian 
net  (see  Appendix  I).  Second,  the  measurements  were  made  using  both  direct 
manometer  board  pressure  readings  and  pressure  transducer  readings  which 
were  sequentially  converted  to  analog  voltages  and  displayed  on  a  digital 
voltmeter.  Third,  at  each  mean  angle  of  attack  setting  made  during  the 
oscillatory  tests,  the  thermally  soaked  differential  pressure  transducer 
outputs  were  also  used  to  check  previously  obtained  results.  The  steady- 
state  normal  force  and  pitching  moment  coefficients  versus  angle  of  attack 
shown  in  Figures  18  and  19  are  the  final  result  of  all  these  measurements. 
The  fact  that  the  normal  force  at  a  =  0  degrees  was  not  zero  for  the 
symmetrical  NACA  0012  airfoil  merely  reflects  the  fact  that  the  mean  air¬ 
flow  angle  in  the  tunnel  differed  from  the  local  horizontal  by  the  noted 
angle  of  zero  lift  (aoj.  =  0.6  degrees).  All  angles  used  in  this  report 
are  referenced  to  local  horizontal,  not  the  angle  of  zero  lift. 

NORMAL  FORCE  AND  MOMENT  HYSTERESIS  LOOPS 


The  smoothed  data  described  in  an  earlier  section  (e.g.,  that  shown 
in  Figure  16)  were  time  histories  of  the  unsteady  normal  force,  pitching 
moment,  and  angle  of  attack.  Additional  information  can  be  obtained  from 
such  data  when  explicit  time  dependency  is  removed.  This  was  accomplished 
conveniently  by  replotting  the  normal  force  and  moment  as  functions  of 
instantaneous  angle  of  attack.  Since  one  complet  cycle  of  data  is  con¬ 
sidered,  the  angle -of -attack  variation  begins  and  ends  at  the  same  point, 
and  the  resulting  curves  are  closed.  The  nature  of  these  closed  curves  is 
best  understood  by  considering  the  case  of  sinusoidal  motion  first. 
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Sinusoidal  Time  Histories 

When  an  airfoil  is  sinusoidally  oscillated  quasi -steadily  over  a 
small  amplitude  range  at  zero  mean  angle  of  attack,  the  resultant  normal 
force  and  moment  versus  angle  of  attack  retrace  the  steady-state  force  and 
moment  curves  for  that  airfoil.  Typical  data  of  this  sort  are  shown  in 
Figure  20a.  If  only  the  frequency  is  increased,  a  phase  shift  between  the 
force  or  moment  and  the  angle  of  attack  begins  to  manifest  itself.  Now 
the  unsteady  force  and  moment  versus  angle  of  attack  curves  enclose  a 
finite  area  as  depicted  in  Figure  20b  (arrows  denote  the  direction  of 
increasing  time).  When  the  mean  angle  of  attack  of  the  motion  is  increased 
such  that  the  maximum  angle  of  attack  of  the  oscillation  exceeds  the 
steady-state  stall  angle,  the  stall  phenomenon  is  delayed.  This  behavior 
is  shown  in  Figure  20c,  where  no  stall  actually  occurs  despite  the  fact 
that  the  steady-state  stall  angle  is  exceeded.  A  further  increase  in  mean 
angle  of  attack  results  in  a  deep  stall  and  recovery  usually  referred  to 
as  dynamic  stall.  Figure  20d  shows  force  and  moment  loops  of  this  latter 
kind.  In  the  classical  or  potential  flow  region  the  unsteady  loops  are 
elliptical,  while  in  the  separated  flow  region  they  can  be  severely 
distorted.  This  discussion  represents  an  oversimplification  of  the 
dynamic  stall  phenomenon.  Actually  the  magnitude  and  character  of  the 
dynamic  stall  delay  are  dependent  on  the  reduced  frequency,  the  penetra¬ 
tion  beyond  static  stall,  and  other  related  factors.  For  a  more  complete 
understanding,  the  reader  is  referred  to  the  previous  literature  on  this 
subject  (References  7,  8,  9,  11,  18,  19,  20,  and  21,  to  name  but  a  few). 

Additional  knowledge  of  the  behavior  of  oscillating  airfoils  can  be 
gained  from  observing  the  manner  in  which  moment  loop  areas  are  enclosed. 
Reference  20  has  shown  that  aerodynamic  damping  of  an  airfoil  is  propor¬ 
tional  to  the  total  area  enclosed  by  these  moment  loops.  Specifically, 
the  following  expression  for  damping  is  taken  from  this  reference: 


M 
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two-dimensional  aerodynamic  damping 
area  enclosed  by  the  moment  loop  and 
oscillatory  amplitude 

It  was  also  shown  in  Reference  20  that  counterclockwise  enclosures  of  the 
moment  loop  are  stabilizing  while  clockwise  enclosures  are  destabilizing. 
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Thus,  as  the  mean  angle  of  attack  is  increased  from  Figure  20c  to  20d,  it 
can  be  seen  that  a  portion  of  the  moment  loop  becomes  unstable  (i.e.,  the 
lower  portion  of  the  figure  eight  loop  in  Figure  20d  is  clockwise)  and  the 
overall  aerodynamic  deunping  is  decreased.  The  aerodynamic  damping  is  of 
the  utmost  importance  in  investigations  of  dynamic  stall  or  stall  flutter. 

The  characteristics  of  the  normal  force  and  moment  hysteresis  loops 
discussed  above  were  found  consistently  throughout  the  sinusoidal  tests 
conducted  in  this  program. 

Nonsinusoidal  Time  Histories 

The  cam-follower  system  used  in  these  tests  made  it  possible  to  obtain 
a  wide  variety  of  airfoil  angle -of -attack  time  histories.  An  extreme 
example  of  nonsinusoidal  motion  was  obtained  from  a  ramp  cam,  which  pro¬ 
duced  regions  of  constant  or  nearly  constant  angular  velocity  connected  by 
small  regions  in  which  the  direction  of  motion  was  reversed.  Figure  21  is 
a  plot  of  the  theoretical  nondimensional  time  history  of  the  angle  of 
attack  and  its  first  two  time  derivatives  as  designed  into  the  ramp  cam. 

The  solid  curves  represent  the  forward  motion  of  the  cam.  Under  these 
circumstances  the  angular  velocity  during  the  upstroke  was  twice  the 
angular  velocity  during  the  downstroke.  This  cam  could  also  be  operated 
in  reverse  (denoted  by  backward  ramp  motion),  which  is  depicted  in 
Figure  21  by  the  dashed  curves.  Here  the  angular  velocity  during  the 
downstroke  was  twice  the  angular  velocity  during  the  upstroke.  The  angular 
acceleration  was  zero  everywhere  except  over  the  regions  of  maximum  and 
minimum  angle  of  attack,  where  it  pulsed  rapidly  to  a  minimum  and  a  maximum 
value,  respectively. 

One  objective  of  this  portion  of  the  experiment  was  to  examine  the 
effect  of  varying  the  time  rate  of  change  of  a  from  one  part  of  the  cycle 
to  another.  This  is  indicated  in  Table  VI,  in  which  values  of  A  are 
listed  for  the  constant  portions  of  the  motions  for  the  four  frequencies 
tested,  for  both  forward  and  backward  motions.  For  example,  it  is  seen 
that  the  forward  cam  at  10  cps  and  the  backward  cam  at  20  cps  will  have 
identical  values  of  A  over  the  upstroke  but  will  have  different  values 
over  the  downstroke.  The  normal  force  response  for  these  conditions  is 
shown  in  Figure  22.  The  solid  curve  labeled  F-10  was  measured  with  the 
forward  cam  at  10  cps,  and  the  dashed  curve  labeled  B-20  was  measured  with 
the  backward  cam  at  20  cps,  for  oscillations  of  *8  degrees  about  a  mean 
angle  of  attack  of  11  degrees.  It  is  seen  that  during  the  upstroke,  both 
of  these  loops  coincide  over  most  of  the  a  range.  However,  afte.’  the 
maximum  angle  of  attack  was  reached,  the  angular  velocity  on  the  downstroke 


was  four  times  greater  for  the  backward  ramp  (B-20)  than  for  the  forward 
ramp  (F-10).  The  higher  negative  A  in  the  backward  ramp  appears  to  have 
caused  the  flow  to  remain  separated  over  a  wider  range  of  angle  of  attack 
than  the  smaller  negative  A  of  the  forward  ramp,  although  a  small  portion 
of  this  difference  between  the  two  curves  may  be  attributable  to  increased 
phase  shift  due  to  a  larger  value  of  reduced  frequency.  Conversely,  the 
lower  negative  value  of  A  on  the  dcwnstroke  of  the  forward  ramp  appears  to 
have  allowed  the  flow  to  reattach  close  to  the  mean  angle  of  attack. 

A  different  trend  can  be  observed  on  the  upstroke  by  comparing  two 
loops  for  two  different  frequencies  for  the  forward  ramp.  In  addition  to 
the  10-cps  condition  (F-10),  another  solid  curve  has  been  plotted  for  a 
20-cps  oscillation  of  the  forward  ramp  about  a  mean  angle  of  attack  of 
11  degrees  (F-20).  Here  the  doubling  of  A  at  20  cps  (F-20)  relative  to  the 
10-cps  condition  (F-10)  for  both  the  upstroke  and  the  dcwnstroke  caused  an 
overall  enlargement  of  the  loop;  the  maximum  normal  force  attained  was 
greater,  and  the  delay  in  normal  force  recovery  over  the  dcwnstroke  was 
also  greater  for  curve  F-20  relative  to  curve  F-10.  Additional  cases  for 
normal  force  were  examined  and  they  showed  a  comparable  trend,  consistent 
with  the  results  shown  above. 

Similar  families  of  loops  were  obtained  for  pitching  moment.  A  sample 
is  shown  in  Figure  23  for  the  forward  ramp  at  10  cps  (solid  curve,  F-10) 
and  the  backward  ramp  at  20  cps  (dashed  curve,  B-20).  Once  again,  the 
curves  coincide  over  the  region  where  the  values  of  A  are  identical  and 
diverge  where  the  values  of  A  are  not  equal.  However,  the  differences  on 
the  return  stroke  are  not  so  accentuated  here  as  in  the  case  of  normal 
force,  until  after  recovery  from  stall  appears  to  have  taken  place. 

Similar  trends  were  observed  for  the  moment  curves  obtained  at  other 
frequencies. 

An  additional  nonsinusoidal  cam  was  built  to  simulate  the  cyclic 
pitch  variation  of  a  helicopter  rotor  blade.  However,  when  tests  were  run 
with  this  cam,  it  was  found  that  surface  irregularities  caused  the  airfoil 
to  respond  erratically  and  imposed  large  loads  on  the  follower  bearings, 
producing  results  which  were  difficult  to  interpret.  These  tests  were 
discontinued  after  a  short  time ,  and  other  tests  with  the  sinusoidal  cam 
were  substituted. 

UNSTEADY  ffiESSUBE  INVESTIGATION 

During  the  rigid  sinusoidal  portion  of  the  test  program,  a  few  cases 
were  run  at  a  mean  angle  of  attack  of  11  degrees  and  at  frequencies  of 
50  and  75  cps  in  which  single  surface  pressure  measurements  were  made 
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(see  item  4  in  Table  II ).  These  data  were  obtained  by  taping  the  upper 
and  lower  surfaces  of  the  airfoil,  respectively,  so  that  the  differential 
pressure  transducers  were  sensing  the  individual  surface  pressures.  The 
purpose  of  this  test  was  to  establish  whether  a  wave  type  fluctuation  in 
pressure  could  be  observed  on  either  or  both  sides  of  the  airfoil  during 
the  oscillations. 

Figure  24  shows  the  original  chordwise  variation  in  differential 
pressure  measured  directly  from  the  transducers  with  both  sides  active 
compared  with  the  same  differential  distribution  reconstructed  from  the 
single  surface  measurements,  both  taken  at  an  instantaneous  incidence 
angle  of  3  degrees  during  a  sinusoidal  oscillation.  Note  that  these 
pressure  differentials  are  measured  relative  to  the  pressure  at  the  mean 
angle  of  attack;  i.e.,  the  Ap  (x)  plotted  here  is  taken  to  be  zero  along 
the  chord  at  «  =  This  result  is  presented  here  to  verify  the  ability 

of  the  single  surface  measurements  to  reproduce  faithfully  the  results 
obtained  from  the  differential  transducers. 

The  pressure  and  suction  surfaces  of  the  airfoil  are  considered 
separately  in  Figures  25  and  26,  respectively.  Here  the  instantaneous 
chordwise  pressure  distributions  on  each  side  of  the  airfoil  are  plotted 
at  various  instantaneous  angles  of  attack  during  a  cycle  of  motion. 

Again  the  pressures  are  measured  relative  to  the  pressure  at  the  mean 
angle  of  attack.  The  solid  lines  represent  the  conditions  of  incidence 
angle  increasing  relative  to  =  11  degrees  (  at  >0),  and  the  dashed  lines 
represent  the  conditions  of  incidence  angle  decreasing  relative  to  = 

11  degrees  (a  <  0) . 

In  Figure  25,  it  is  seen  that  the  pressure  distribution  along  the 
pressure  surface  changed  very  little  over  much  of  the  cycle,  either  in 
magnitude  or  in  shape,  although  it  changed  sign  relative  to  the  mean 
pressure,  as  expected.  In  contrast,  the  pressure  distribution  over  the 
suction  surface  changed  drastically  as  incidence  angle  was  changed  (Figure 
26).  Here  it  is  seen  that  for  the  angles  decreasing  from  the  mean  angle, 
a  potential  flow  condition  was  established  and  the  pressure  distributions 
are  well  behaved  and  exhibit  the  usual  leading -edge  peak.  However,  for 
the  angles  increasing  from  the  mean  angle,  it  appears  that  a  dynamic  stall 
condition  was  encountered.  Note  that  in  Figure  26  the  potential  behavior 
of  the  curve  for  a  =  L3.7  degrees  becomes  distorted  for  a  =  l6.3  degrees 
and  appears  to  have  developed  a  propagating  wave  for  a  =  19.0  degrees. 
Additional  pressure  curves  were  examined  for  a  decreasing  from  the  peek 
value  of  19  degrees  (not  shown  here  to  avoid  the  confusion  of  overlapping 
curves),  and  it  was  found  that  a  distinct  pressure  wave  could  be  discerned 


propagating  rearward  along  the  chord  until  it  be  cane  too  small  to  identify 
at  approximately  the  40  or  50  percent  chord  station.  Almost  identical 
results  were  observed  for  the  higher  frequency  (f  =  75  cps)  case. 

The  rearward  propagation  of  a  pressure  wave  strongly  suggests  a 
vortex  shedding  phenomenon  as  described  by  Ham  in  Reference  19*  Further¬ 
more,  the  relative  lack  of  response  on  the  pressure  surface  (Figure  25) 
indicates  that  the  bulk  of  the  dynamic  stall  activity  occurs  on  the  suction 
surface.  This  phenomenon  can  also  be  inferred  from  the  figures  in 
Reference  22,  and  it  was  clearly  evident  in  the  motion  picture  (Reference 
23)  of  the  animated  pressure  distributions  which  accompanied  Reference  22. 
Thus,  it  should  be  possible  to  use  even  the  unsteady  differential  pressure 
data  to  examine  some  of  the  details  of  the  dynamic  stall  phenomenon.  This 
work  will  be  undertaken  in  the  future. 


UNSTEADY  AIRFOIL  DATA 


GENERALIZATION  OF  UNSTFADY  SINUSOIDAL  DATA 

The  process  used  to  generalize  these  results  has  already  been  briefly 
described  in  a  previous  section  and  has  been  thoroughly  discussed  in 
Reference  8.  As  in  the  past,  the  normal  force  and  moment  coefficient  data 
were  tabulated  as  functions  of  o,  A,  and  B,  and  these  tabulations  were 
stored  in  the  computer.  (They  are  also  found  in  this  report  in  Tables  IV 
and  V  for  Cn  and  Cm,  respectively.)  The  differences  between  this  tabula¬ 
tion  and  our  previous  tabulation  (Reference  8)  were  many  and  varied.  We 
shall  be  concerned  here  only  with  a  few  of  the  more  important  differences. 
First,  the  model  scale  wa .«  considerably  different;  this  time  the  chord  was 
5  inches  and  in  Reference  8  it  was  24  inches.  Another  important  difference 
lay  in  the  method  of  extracting  points  from  the  original  time  histories. 
Previously,  a  single  cycle,  chosen  at  random,  was  analyzed  and  manual 
procedures  were  used  to  interpolate  the  normal  force  and  moment  values, 
cross-plot  the  results,  and  extract  the  final  values  for  tabulation.  In 
the  present  study,  the  single  cycle  was  replaced  by  a  time  -  averaged  cycle 
and  most  of  the  intermediate  steps  were  performed  by  a  computer.  This  led 
to  a  self-consistent  set  of  data  which  was  smoother  than  that  obtained 
previously.  Finally,  the  parameters  chosen  for  this  experiment  were 
specifically  picked  to  extc-rd  the  useful  tabulated  range  beyond  that  of 
the  previous  table. 

It  is  instructive  to  examine  the  a,  A,  B  table  generated  by  the 
present  test  and  to  compare  its  range  with  the  required  parameter  ranges 
of  typical  conditions  experienced  in  flight.  To  accomplish  this,  we  can 
combine  Equations  (12)  and  (13)  of  Reference  8  to  eliminate  the  angular 
displacement  from  the  mean,  and  after  some  manipulations,  we  obtain 
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This  equation  shows  that  contours  of  constant  k  and  a  will  be  ellipses  in 
the  A,  B  plane.  Furthermore,  the  major  and  minor  axes  of  any  given  ellipse 
will  be  directly  governed  by  the  size  of  both  the  reduced  frequency  and 
the  angular  amplitude,  being  large  when  ka  is  large  and  small  when  kft  is 
small.  This  is  clearly  illustrated  in  Figure  27.  The  solid  outside 
ellipse  is  the  limit  of  the  conditions  tested  in  this  study,  and  the  data 


tabulation  (presented  in  Tables  IV  and  V)  is  arranged  roughly  in  the  same 
pattern.  (Note  that  the  moment  data  in  Table  V  extend  farther  into  the 
corners  of  the  cable  than  do  the  normal  force  data  in  Table  IV,  because 
the  Cm  data  were  extrapolated  through  a  series  of  cross -plots  to  increase 
their  useful  range.)  For  comparison,  the  large  solid  rectangle  represents 
the  tabular  limits  of  the  old  data  used  in  Reference  8.  The  dashed  curves 
labeled  a  through  d  are  the  anticipated  limits  for  various  application 
requirements.  Curves  a  and  b  represent  the  bounds  of  typical  helicopter 
rotor  blades  executing  one-per-rev  motions  at  *8  degrees  amplitudes  for 
high  and  low  aspect  ratios,  respectively  (the  differences  in  the  curves 
being  associated  with  differences  in  blade  chord  length).  Curve  c  repre¬ 
sents  a  sinusoidal  stall  flutter  limit  of  eight -per-rev  at  ±4  degrees,  and 
curve  d  represents  the  maximum  frequency  boundary  of  the  nonsinusoidal 
ramp  motion  during  the  present  program.  It  is  clear  from  this  figure  that 
the  sinusoidal  data  obtained  is  adequate  for  the  analysis  of  the  typical 
cases  shown  here . 

one  of  our  major  objectives  in  this  program  was  to  examine  the 
ability  of  the  generalized  data  taken  from  a  sinusoidal  test  to  predict 
the  aerodynamic  response  of  the  airfoil  as  it  executed  a  more  general 
periodic  (nonsinusoidal)  motion.  Therefore,  an  obvious  first  requirement 
of  the  tabulated  data  is  that  for  a  sinusoidal  motion,  the  predicted 
normal  fore*  and  moment  loops  must  faithfully  reproduce  the  original  loops. 
It  is  seen  from  Figures  28  and  29  that  this  condition  is  essentially  sat¬ 
isfied.  These  figures  represent  a  sampling  of  the  normal  force  and 
moment  comparisons  (respectively)  made  between  the  original  data  and  the 
reproduced  loops.  The  left  panels  in  each  figure  show  the  effect  of  mean 
angle -of -attack  variation  for  a  frequency  of  approximately  50  cps.  The 
right  panels  in  each  figure  show  the  effect  of  frequency  at  a  constant 
mean  angle  of  attack  of  11  degrees.  In  almost  all  cases,  the  agreement 
between  the  original  data  and  the  reproduced  loops  is  excellent, 
particularly  at  the  higher  frequencies.  In  a  few  instances,  the  agreement 
at  low  frequency  shows  some  deterioration,  but  it  is  still  acceptable  for 
our  present  purposes.  This  is  shown  in  Figure  29c,  which  is  typical  of 
the  maximum  deviation  of  a  reproduced  loop  from  the  original  observed 
during  this  study.  Thus,  the  present  force  and  moment  tabulation  has  been 
shown  to  be  an  extremely  accurate  duplication  of  the  experimental  data. 

This  accuracy  appears  to  be  significantly  greater  than  that  demonstrated 
by  any  other  method  for  which  results  have  been  presented  in  the 
literature . 
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A  few  cases  of  sinusoidal' motion  were  also  run  at  an  oscillatory 
amplitude  of  *3  degrees.  Once  again,  the  unsteady  normal  force  and  moment 
were  predicted  using  the  a,  A,  and  B  tabulation  generated  from  the  -8 
degree  sinusoidal  tests.  Figures  30  and  31  show  the  results  of  comparing 
the  original  -3  degree  experimental  normal  force  and  moment  to  that 
predicted  by  the  tabulation.  Both  figures  are  for  a  mean  angle  of  attack 
of  14  degrees.  Figure  30  is  for  an  oscillatory  frequency  of  12.6  cps  and 
shows  very  good  agreement  in  both  normal  force  and  pitching  moment. 

Figure  31  is  for  an  oscillatory  frequency  of  31.1  cps,,  and  the  agreement 

is  only  fair  for  the  decreasing  incidence  branch,  although  it  is  good  for 

increasing  incidence. 

USE  OF  DATA  TABULATION  TO  PREDICT  GENERAL  RESPONSE  ' 

It  was  stated  in  an  earlier  section  that  one  object  of  the  ramp 
motion  tests  was  to  examine  the  effect  of  variable  a  within  a  cycle  of 

motion  on  the  blade  aerodynamic  response.  A  more  important  object  was  to 

test  the  validity  of  the  hypothesis  that  sinusoidal  data  could  be  general¬ 
ized  by  the  UAC  a,  A,  B  method  to  predict  the  aerodynamic  .response  of  an 
airfoil  to  a  nonsinusoidal  motion.  The  abrupt  changes  in  angular  velocity 
and  acceleration  provided  an  extreme  environment  in  which  to  test  this 
hypothesis.  ,  , 

I 

Figure  32  is  a  plot  of  the  actual  nondimensional  time  history  of  the 
angle  of  attack  produced  by  the  ramp  cam  when  it  was  rotated  in  the  for¬ 
ward  direction.  When  this  time  history  was  analyzed  for  use  with  the 
tabulation,  it  was  found  that  there  were  minor  deviations  between  the  true 
(or  measured)  angular  ppsition  and  the  ideal  (constant  a  )  trajectory. 
These  irregularities  resulted  from  rmall  errors  in  the  cam  surfaces. 
Although  these  anomalies  were  within  acceptable  limits,  they  did  cause  . 
some  spurious i oscillations  when  the  angle  of  attack  was  numerically 
differentiated  to  obtain  the  angular  rate  term,  A,  and  the  angular  ^ccel,- 
eration  term,  B.  This  problem  wds  avoided  in  the  present  application  by 
fitting  the  angle -of -attack  time  histories  with  high-order  polynomials 
(up  to  fourteenth  ordpr),  which  were  then  differentiated  exactly  rather  ' 
than  numerically.  The  polynomial  fits  produced  the  actual  time  histories 
with  good  accuracy. 

•  l 
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Figures  33  and  3^  show  some  typical  comparisons  between  measured 
force  and  'moment  loops  for  the  forward  ramp  motion  and  those  predicted 
from  the  a,  A,  and  B  tabulation  of  sinusoidal  data.  In  all  cases,  the 
predicted  Cnj^^  was  within  5  percent  of  the  measured  value  in  magnitude 
but  did  tend  to  occur  at  a  slightly  lower  arjgle  of  attaok.  After  stall, 
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the  predicted  normal  force  decreased  more  rapidly  and  tended  to  remain 
stalled  over  a  greater  angle  of  attack  range  than  did  the  measured  normal 
force.  Tlie  overall  agreement  of  normal  force  is  considered  to  be  very 
good.  The  moment  loops  of  Figure  33  also  appear  tc  be  in  gr_j  agreement 
with  regard  to  shape  (including  crossover  points),  although  a  small  shift 
in  magnitude  is  apparent.  Since  stall  flutter  is  primarily  a  torsional 
problem,  the  moment  loops  are  somewhat  more  important  than  the  normal 
force  loops.  A  measure  of  how  well  the  measured  and  predicted  moment  loops 
agree  can  be  obtained  by  comparing  the  aerodynamic  damping  calculated  from 
each  set  of  loops.  An  analysis  of  aerodynamic  damping  is  presented  in 
Reference  20  and  is  discussed  briefly  in  a  preceding  section  of  this 
report.  The  expression  for  aerodynamic  damping  given  by  Equation  (6)  was 
used  to  obtain  the  damping  for  both  the  measured  and  predicted  moment 
loops  for  all  tests  with  the  ramp  time  history  of  Figure  32.  These  results 
are  shown  in  Figure  35*  The  agreement  at  low  reduced  frequency  and  low 
mean  angle  of  attack  is  excellent.  At  higher  reduced  frequencies,  the 
damping  for  the  predicted  moment  loops  is  consistently  less  (i.e.,  more 
negative  or  tending  more  toward  negative)  than  the  damping  for  measured 
moment  loops,  particularly  for  oscillatory  motions  which  penetrate  and 
return  from  stalled  flow.  Similar  results  were  obtained  for  the  backward 
ramp,  as  shown  in  Figure  36  for  the  normal  force,  in  Figure  37  for  the 
moment,  and  in  Figure  38  for  the  aerodynamic  damping. 

FLEXURED  TEST  RESULTS 


The  bulk  of  the  tests  conducted  in  this  program  were  with  a  rigid 
coupler  between  the  transmission  shaft  and  the  airfoil.  It  was  possible 
to  simulate  the  stall  flutter  problem  two -dimensionally  by  replacing  this 
rigid  coupler  with  flexible  ones.  Three  different  Bendix  flexures  were 
used  which  resulted  in  natural  frequencies  of  the  airfoil  of  41.5,  52  and 
75  cps.  It  was  stated  earlier  that  the  airfoil  was  constructed  in  two 
parts  and  was  mounted  in  thrust  bearings  in  the  sidewalls  of  the  wind 
tunnel.  These  thrust  bearings,  which  were  required  to  hold  the  two  model 
parts  together,  introduced  more  mechanical  damping  in  the  system  than  is 
normally  desirable  for  a  stall  flutter  test.  Typical  values  of  the  no -flew 
damping  obtained  from  the  decay  rate  of  the  plucked  spring -mounted  airfoil 
were  between  3.9  and  10.2  percent  of  critical  damping.  Despite  this 
amount  of  damping,  it  was  felt  that  a  stall  flutter  test  was  still  meaning¬ 
ful,  particularly  when  it  was  observed  that  significant  flutter  response 
still  occurred. 


The  airfoil  as  described  was  oscillated  sinusoidally  at  an  amplitude 
of  -8  degrees  about  mean  angles  of  attack  of  6,  11,  and  16  degrees.  At 
each  mean  angle  of  attack,  three  different  excitation  drive  frequencies 
were  used  to  simulate  the  one -per -rev  cyclic  pitch  of  a  helicopter  blade. 
The  combinations  of  natural  frequencies  and  excitation  drive  frequencies 
resulted  in  a  range  of  frequency  ratios  from  2.76  to  10.0  cycles -per -rev, 
which  is  within  the  typical  range  of  most  helicopters.  The  oscillatory 
amplitude  of  *8  degrees  was  also  representative  of  typical  helicopter 
motion . 

Use  of  flexures  did  not  change  any  of  the  test  procedures  described 
earlier.  Both  position  gages  were  required  for  these  tests;  one  measured 
the  excitation  output  of  the  transmission  ( one -per -rev ) ,  while  the  other 
measured  the  response  of  the  airfoil  (many -per -rev) . 

In  addition  to  the  direct  measurement  of  the  airfoil  response,  an 
analysis  was  performed  in  which  the  generalized  a,  A,  B  tabulation  for  the 
moment  (Table  V)  was  used  to  predict  the  airfoil  response.  Details  of  the 
numerical  procedure  used  are  given  in  Appendix  III.  The  various  physical 
characteristics  of  the  airfoil  required  in  this  analysis  are  given  in 
Table  VTI.  Several  values  of  the  logarithmic  decrement  of  the  mechanical 
damping  were  used  in  the  analysis  since  the  actual  value  of  this  damping 
during  testing  was  not  known.  This  was  because  the  differential  pressure 
across  the  sidewalls  during  the  tests  caused  additional  thrust  loadings 
on  the  bearings,  which  in  turn  caused  an  unknown  increase  in  mechanical 
damping.  However,  it  was  felt  that  the  total  mechanical  damping  of  the 
system  did  not  exceed  the  values  used  in  the  analysis. 

Some  typical  results  of  these  tests  are  shown  in  Figures  39  and  1+0 
for  different  combinations  of  frequency  and  angle  of  attack.  The  solid 
curves  denote  the  experimental  data,  while  the  dashed  curves  denote  the 
results  of  the  numerical  solution  for  the  same  conditions  using  the 
tabulated  values  of  Cm.  Each  of  the  test  results  shown  is  one  cycle  of 
time -averaged  motion.  The  time -averaging  process  was  applied  to  Uo  cycles 
of  motion  as  in  other  parts  of  this  study.  To  verify  the  fact  that  this 
averaging  procedure  did  not  wash  out  important  features  of  the  motion,  it 
was  also  done  for  five  cycles,  and  no  essential  difference  was  found. 

Although  the  agreement  between  test  and  analysis  is  not  perfect,  the 
results  indicate  that  the  a,  A,  B  method  does  yield  results  which  contain 
all  the  essential  characteristics  of  the  test  results.  Each  graph  in 
Figures  39  and  Uo  shows  three  pairs  of  curves,  one  pair  for  each  mean 
angle  of  attack.  For  convenience,  each  curve  was  plotted  as  a  deviation 
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from  the  mean  angle  of  attack,  and  on  these  figures  each  pair  of  curves  is 
drawn  from  its  own  origin  on  the  ordinate.  Figure  39  plots  three  differ¬ 
ent  excitation  frequencies  for  an  airfoil  natural  frequency  of  41.5  cps, 
while  Figure  40  shows  the  same  three  excitation  frequencies  for  an  airfoil 
natural  frequency  of  52  cps.  In  all  six  plots,  it  should  be  observed  that 
the  occurrences  of  natural  blade  oscillations  at  high  angle  of  attack 
which  typify  stall  flutter  were  all  predicted.  Note  also  that  the  initial 
cycle  of  predicted  blade  response  always  coincided  with  the  measured  res¬ 
ponse.  Since  these  blade  oscillations  and  their  inception  points  are  the 
most  important  aspects  of  stall  flutter,  these  results  are  particularly 
significant.  At  the  lowest  angle  of  attack  (cc^  =  6  deg)  and  the  lower 
blade  natural  frequency  (41.5  cps)  shown  in  Figure  39>  the  analytical  sol¬ 
ution  behaved  almost  sinusoidally  while  the  test  data  exhibited  a  broadened 
cycle  on  the  upstroke  (i.e.,  the  test  angle -of -attack  time  history  on  the 
upper  half  of  the  cycle  extended  over  a  time  greater  than  one -half  the 
period  of  the  excitation  frequency).  This  phenomenon  was  not  nearly  as 
severe  for  the  aM  =  6  degrees  cases  at  an  airfoil  natural  frequency  of 
52  cps  (see  Figure  4o).  For  all  other  cases,  the  agreement  between  theory 
and  test  results  was  considered  to  be  excellent. 
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PROCESSING  OF  THE  UNSTEADY  DATA  FOR  USE  IN 
HELICOPTER  CONTROL  LOAD  CORRELATION  STUDIES 


MODIFICATION  OF  THE  ORIGINAL  UNSTEADY  DATA  TABLES 


The  unsteady  data  tabulations  presented  in  Reference  8  were 
reexamined,  and  modifications  were  made  in  an  attempt  to  improve  correla¬ 
tion  with  hysteresis  loops  for  sinusoidal  motion.  The  procedure  used  was 
to  examine  the  cross -plots  of  normal  force  and  pitching  moment  coefficients 
versus  a,  A,  ana  B  and  to  smooth  out  irregular  variations.  Examples  of  the 
original  and  the  smoothed  cross -plots  are  shown  in  Figures  4l  and  42.  The 
effect  of  these  changes  on  the  ability  to  reconstruct  the  original 
hysteresis  loops  was  evaluated,  and  those  modifications  which  led  to 
general  improvement  in  correlation  were  incorporated  into  the  tables. 
Because  of  the  scatter  in  the  data,  alterations  which  improved  correlation 
of  some  loops  at  times  deteriorated  that  of  others .  Samples  of  some  of 
the  improvements  in  correlation  are  shown  in  Figure  43. 

SCALING  OF  THE  NEW  UNSTEADY  DATA 


Application  of  the  unsteady  data  to  rotor  blade  operating  conditions 
requires  that  the  measured  data  for  Mach  number  0.325  be  scaled  to  obiain 
approximations  of  higher  Mach  number  unsteady  characteristics.  As  in 
Reference  8,  scaling  procedures  based  on  the  steady-state  stall  angle  were 
chosen  which  transformed  the  steady-state  Cn  and  Cm  data  obtained  for  the 
wind  tunnel  model  to  steady  normal  force  and  moment  curves  for  a  0012  air¬ 
foil  over  a  range  of  Mach  numbers.  These  steady-state  scaling  laws  were 
then  assumed  to  apply  to  the  unsteady  blade  characteristics  as  well.  It 
was  shown  in  Reference  8  that  normal  force  and  pitching  moment  coefficient 
hysteresis  loops  measured  at  M  =  0.4  and  M  =  0.6  could  be  predicted  with 
reasonable  accuracy  by  scaling  the  measured  data  for  M  =  0.2.  The  agree¬ 
ment  between  measured  and  scaled  data  was  good  for  M  =  0.4  but  was  less 
favorable  for  M  =  0.6.  Although  it  would  be  desirable  to  have  experimental 
unsteady  data  at  high  Mach  numbers,  it  is  felt  that  the  present  data  can 
be  adequately  scaled  to  the  Mach  numbers  typical  of  the  retreating -blade 
portion  of  the  rotor  disk  where  stall -related  unsteady  effects  are  signif¬ 
icant. 

Other  factors  included  in  the  data  transformation  process  are  a 
correction  for  the  finite  effective  aspect  ratio  of  the  model  and  a 
correction  for  the  slight  difference  between  the  chosen  oc  =  0  degrees 
reference  line  and  the  measured  angle  of  zero  normal  force.  For  the 
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steady-state  model  data  shown  In  Figure  18,  the  normal  force  curve  slope 
is  only  about  0. 0836/degree  and  the  measured  angle  of  zero  normal  force 
is  0.6  degrees.  The  treatment  of  these  two  factors  led  to  a  data  scaling 
technique  slightly  different  from  that  used  in  Reference  8.  This  will  now 
be  discussed. 

Normal  Force 

To  scale  the  measured  unsteady  normal  force  data  for  Mach  number 
effects  using  a  process  based  on  steady-state  stall  angle,  it  was  first 
necessary  to  replot  the  steady  model  Cn  data  versus  the  relative  stall 
angle  parameter 
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The  normal  force  stall  angle  or  point  of  departure  from  a  linear  normal 
force  curve,  asn,  and  the  angle  of  zero  normal  force,  «ol>  are  taken  from 
Figure  18  to  be  12.0  degrees  and  0.6  degree,  respectively.  The  resulting 
plot,  shown  in  Figure  44,  is  then  the  basic  curve  from  which  all  normal 
force  data  are  scaled.  The  unsteady  data  tables  used  in  the  following 
correlation  studies  also  use  this  calculated  <rn  as  an  independent  variable 
in  place  of  measured  angle  of  attack. 

The  process  for  mapping  the  steady  model  data  onto  an  arbitrary 
steady  normal  force  curve  is  depicted  in  Figure  45  and  is  described  as 
follows  : 

(1)  Choose  the  point  of  departure  from  linearity,  asn,  and  re plot 
the  given  normal  force  curve  versus  <xn  =  a/asn- 

(2)  Calculate  the  ratio  Kn  where  Kn  =  Cnstall/Cnstall  model 
where  Cngtall  =  Cn(asn). 

(3)  Plot  an  intermediate  scaling  curve  Cn'(<rn)  where  ^n'^n)  = 

^  ^nmode 1 (°n ) • 

(4)  Determine  an  additive  correction  curve  8n(<rn)  such  that 
cn(°h)  =  cn '  (°n )  +  8n(<rn). 
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This  procedure  was  carried  out  for  full-scale  NACA  0012  Cn  data.  The 
steady-state  data  taken  from  Reference  24  were  obtained  from  a  Sikorsky- 
built  16 -inch-chord  airfoil  section  which  was  tested  in  the  UAC  8-foot 
Subsonic  Wind  Tunnel  at  Reynolds  numbers  from  2.6  to  5.3  million.  The  cn 
characteristics  are  shown  in  Figure  46.  The  scaling  curves  used  to  trans¬ 
form  the  model  data  to  the  full-scale  curves  are  shown  in  Figure  47. 

After  the  three  quantities  agn,  Kn  and  8n(crn)  are  chosen  for  each 
Mach  number  and  the  assumption  is  made  that  the  same  scaling  laws  can  be 
applied  to  the  unsteady  normal  force  coefficients,  the  following  computa¬ 
tion  procedure  is  used  to  arrive  at  unsteady  Cn  coefficients : 

(1)  Compute  the  rotor  blade  section  conditions  <rn,  M,  A,  B. 

(2)  Compute  the  dimensionless  stall  angle  parameter,  trn  =  a/asn* 

(3)  Enter  the  tabulation  of  the  unsteady  Cn  coefficients  with  <rn, 

A,  B  and  extract  Cnmodel(<rn,  A,  B). 

(4)  Calculate  Cn'(<rn,  A,  B)  =  Kn(M)  A,  B). 

(5)  Calculate  Cn(crn,  A,  B)  =  Cn‘(<rn,  A,  B)  +  8n(<rn). 

The  scaling  parameters  used  in  this  Mach  number  scaling  procedure  are 
given  in  Table  VIII. 

Pitching  Moment 

As  in  the  case  of  the  normal  force,  the  Mach  number  scaling  procedure 
used  for  the  pitching  moment  is  based  on  the  steady-state  stall  angle. 
Figure  48  shows  the  NACA  0012  steady-state  pitching  moment  curves  used  in 
Reference  8.  For  each  Mach  number,  the  pseudo-stall  angle  parameter,  arsm, 
is  chosen  to  scale  the  measured  M  =  0.325  model  data  to  the  appropriate 
full-scale  curve.  The  dimensionless  pitching  moment  parameter 
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defines  the  relative  degree  of  pitching  moment  stall  of  the  model.  From 
Figure  19,  asm  is  taken  to  be  l4  degrees,  while  =  0.6  degree  as  before, 
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In  a  manner  analogous  to  that  used  for  the  normal  force  coefficients, 
the  parameter  asm  and  an  additive  pitching  moment  correction  curve  8m(o-m) 
were  chosen  to  scale  the  full-scale  Cm  characteristics  from  the  steady 
model  data.  The  best  values  of  asm  and  8m(<rm)  for  each  Mach  number  are 
chosen  as  follows : 

(1)  Choose  the  stall  angle  parameter,  asm,  which  best  maps  the  stall 
region  of  the  full-scale  curve  from  the  model  data,  and  replot 
the  pitching  moment  curve  versus  crm  =  «/asm. 

(2)  Determine  the  additive  correction  curve  8ra(crm)  such  that 
cm  (°m )  =  cmm0del^am^  +8m(°in)* 

Figure  49  shows  some  of  the  results  obtained  with  this  procedure. 

The  asm  and  8m(<rm)  values  were  used  in  the  following  calculation  procedure 
for  pitching  moment  coefficient : 

(1)  Compute  the  following  two-dimensional  conditions  at  a  rotor 
blade  section  a,  M,  A,  B. 

(2)  Compute  the  dimensionless  stall  angle  parameter  frm  =  a/asm 
where  arsm  corresponds  to  the  blade  section  Mach  number  M. 

(3)  Enter  the  pitching  moment  tabulation  with  trm,  A,  B  and  extract 
Cmmodel  5  A’  B^' 

(4)  Calculate  q„(<rm,  A,  B)  =  A>  B)  +M(rm)' 

Table  VIII  lists  the  moment  stall  angles,  asm,  used  in  the  present 
calculations . 
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CORRELATION  STUDY 


The  purpose  of  the  correlation  study  was  to  evaluate  the  analytical 
capability  for  predicting  stall -flutter -induced  control  loads  using  the 
unsteady  aerodynamic  model  described  herein.  Three  rotor  systems  having 
large  retreating  blade  control  loads  were  chosen  for  this  study:  the 
NH-3A  (S-61F)  compound  helicopter  rotor  at  high  speed,  the  CH-53A  helicop¬ 
ter  rotor  under  heavy  loading,  and  a  dynamically  scaled  model  rotor.  Each 
rotor  system  had  been  extensively  instrumented  to  measure  the  required 
blade  loads  and  responses.  References  25,  26,  and  27  are  the  sources  of 
the  test  data  for  the  NH-3A,  the  CH-53A,  and  the  model  rotor,  respectively. 
The  cases  chosen  for  the  present  study  are  summarized  in  Table  IX  with 
designations  from  the  appropriate  reference. 

The  Normal  Modes  Blade  Aeroelastic  Analysis  was  used  to  calculate  the 
blade  elastic  response  for  each  of  the  flight  conditions.  This  program, 
which  is  described  in  Reference  28,  represents  blade  flatwise,  edgewise, 
and  torsional  elastic  deformation  by  a  summation  of  normal  mode  responses. 

A  timewise  integration  of  the  equations  of  motion  for  a  coupled  rotor 
aerodynamic  and  dynamic  system  is  used  to  solve  for  the  modal  response. 
Using  this  analysis,  the  effect  of  unsteady  aerodynamics  on  the  blade 
dynamic  response  was  studied.  As  detailed  in  Reference  8,  the  unsteady 
airfoil  data  were  incorporated  into  the  blade  response  calculations.  At 
each  azimuthal  and  radial  blade  station  the  blade  angle  of  attack  and  its 
first  two  time  derivatives  were  computed  and  used  to  look  up  the  Cn  and  Cm 
coefficients  from  the  unsteady  tables.  These  coefficients  were  then 
scaled  to  the  local  blade  Mach  number  using  the  procedure  outlined  in  the 
previous  section.  The  importance  of  variable  rotor  inflow  was  evaluated 
using  the  UAC  Prescribed  Wake  Analysis,  which  calculates  the  induced 
velocity  distribution  over  the  rotor  disk  due  to  the  rotor  shed  wake. 

This  procedure  is  described  in  Reference  29.  The  results  of  the  correla¬ 
tion  study  should  indicate  the  importance  of  unsteady  aerodynamic  and 
variable  inflow  effects  in  modelling  rotor  blade  torsional  excitation. 

NH-3A  FLIGHT  TEST  CORRELATION 

The  five-biaded  NH-3A  compound  helicopter  exhibits  pronounced  high- 
frequency  torsional  oscillations  over  a  range  of  high-speed  flight  condi¬ 
tions.  Four  conditions  were  chosen  for  investigation:  a  186-knot  case 
with  C t/o-  =  0.0579>  two  cases  at  approximately  160  knots  and  Crp/cr  =  0.073 > 
and  a  138-knot  60-degree  angle  of  bank  coordinated  turn  at  Crp/tr  =  0.085*+. 
These  cases  were  chosen  because  their  Cx/cr  values  were  among  the  highest 
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experienced  for  this  aircraft.  Because  of  the  auxiliary  propulsion  from 
the  jets,  the  rotor  trim  conditions  are  generally  not  typical  of  pure 
helicopter  operation.  Table  DC  shows  small  tip -path -plane  inclinations 
for  two  of  the  correlation  cases  and  suggests  that  wake  effects  may  have 
an  important  influence  upon  blade  response.  For  a  small  tip -path -plane 
inclination,  the  wake  remains  close  to  the  rotor  and  hence  can  affect  it 
more  than  it  would  for  a  large  tip -path -plane  inclination. 

The  rotor  blades  were  modelled  using  rigid -body  flapping  and  lagging 
degrees  of  freedom  and  five  flatwise,  two  edgewise,  and  two  torsional 
elastic  modes.  In  all  cases,  rotor  shaft  angle  and  cyclic  pitch  were  taken 
from  test  measurements,  and  collective  pitch  was  incremented  to  match 
measured  lift.  The  use  of  two-dimensional  lifting  line  theory  has  con¬ 
sistently  underpredicted  rotor  lift,  making  it  necessary  to  augment 
collective  pitch  in  the  analytical  solution. 

EFFECT  OF  MODIFICATIONS  TO  ORIGINAL  UNSTEADY  DATA 

The  investigation  was  begun  by  running  NH-3A  case  43  with  the  original 
unsteady  aerodynamic  tables  (Reference  8)  which  were  modified  to  smooth 
regions  of  irregular  data.  In  this  calculation,  a  small  fourth -quadrant 
oscillation  was  predicted.  Figure  50  shows  the  analytical  response  with 
and  without  the  data  modifications.  The  level  of  calculated  response  is 
lower  than  the  measured  value  and  does  not  initiate  over  the  nose  of  the 
aircraft . 

The  possibility  that  rotor  trim  might  be  critical  in  determining  stall 
flutter  susceptibility  was  examined.  The  computer  program  was  run  with 
incremental  changes  in  cyclic  and  collective  pitch.  These  cases  generally 
resulted  in  calculated  lifts  and  propulsive  forces  different  from  test 
data.  Figure  51  demonstrates  that  variations  in  longitudinal  cyclic  pitch 
do  modify  the  fourth-quadrant  torsional  response  for  constant  C.p/cr  but  do 
not  change  the  basic  character  of  the  results.  Figure  52  shows  that 
increasing  collective  pitch  can  induce  higher  levels  of  analytic  response 
but  not  the  proper  azimuthal  distribution.  In  this  case,  Gp/cr  varied  from 
0.0475  at  #75  =  7.29  degrees  to  0.067  at  0<j ^  =  9-29  degrees,  as  shown. 

In  the  compound  mode,  nonuniform  rotor  inflow  can  be  expected  to  be  a 
significant  factor  in  determining  blade  angles  of  attack  and,  therefore, 
aerodynamic  forces  and  moments.  Three  forms  of  nonuniform  rotor  inflow 
were  examined.  First,  the  UAC  Prescribed  Wake  Analysis  was  used  to  calcu¬ 
late  the  induced  velocities  caused  by  a  nondistorted  helical  wake.  When 
case  43  was  lerun  with  this  inflow,  there  was  an  improvement  both  in 
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amplitude  and  azimuthal  character  correlation.  As  seen  in  Figure  53  (note 
change  in  scale  in  this  figure),  the  unsteady  aerodynamic  model  predicts  a 
torsional  oscillation  beginning  at  about  azimuth  angle  i jj  =  loO  degrees  and 
persisting  over  the  retreating  side  of  the  rotor  disk.  Second,  the 
Prescribed  Wake  Analysis  was  used  to  calculate  a  rotor  inflow  in  which  the 
distortions  of  the  rotor  wake  due  to  interactions  between  vortex  elements 
were  included.  The  resulting  normal  modes  solution  for  case  43  using  this 
distorted  wake  inflow  was  not  significantly  different  from  the  nondistorted 
wake  solution.  Examination  of  the  small  change  in  calculated  trajectories 
of  tip  filaments  indicates  that  distortion  does  not  significantly  affect 
rotor  inflow  for  this  thrust,  advance  ratio,  and  tip-path -plane  inclina¬ 
tion.  Finally,  the  effect  of  the  aircraft  wing  on  induced  velocity  in  the 
plane  of  the  rotor  was  examined.  Wing  lift  was  estimated  from  aircraft 
gross  weight  and  main  rotor  shaft  tension.  The  UAC  Rotor  Wake  Interference 
Program  was  used  to  describe  induced  velocities  in  the  plane  of  the  rotor 
caused  by  -the  wing  circulation.  The  effect  on  rotor  inflow  was  determined 
to  be  small  for  the  low -wing  S-61F  configuration.  The  torsional  response 
prediction  taking  into  account  distorted  wake  and  wing  inflow  effects  was 
essentially  the  same  as  that  using  the  nondistorted  wake. 

EFFECT  OF  NEW  UNSTEADY  DATA 

When  the  new  unsteady  model  data  became  available,  the  analysis  was 
carried  out  for  NH-3A  cases  28,  39,  43,  and  72.  Figures  54  through  57 
compare  the  steady  and  unsteady  solutions  obtained  for  these  cases  using 
constant  inflow.  The  level  of  torsional  excitation  is  significantly 
underpredicted  by  both  the  steady  and  unsteady  solutions.  The  reason  for 
the  discrepancy  in  the  steady  components  of  the  response  is  not  known. 

This  may  be  the  result  of  changes  in  blade  steady  pitching  moment  intro¬ 
duced  by  trim  tabs  and  tip  weights  used  to  achieve  proper  balancing  and 
tracking.  Also,  there  may  be  a  bias  applied  to  the  data  during  the  sig¬ 
nal  acquisition  and  conditioning  process.  Figure  58  shows  a  sample  of 
flatwise  and  edgewise  stress  correlation  obtained  with  constant  inflow. 

The  level  of  edgewise  vibratory  stress  calculated  with  unsteady  aerody¬ 
namics  and  constant  inflow  is  much  less  than  that  measured  in  test.  Com¬ 
pared  to  calculations  using  steady  aerodynamics,  unsteady  aerodynamic 
solutions  generally  require  lower  collective  pitch  angles  to  match  test 
values  of  lift.  This  results  in  lower  edgewise  drag  loads  and  stresses. 

The  lack  of  reliable  unsteady  drag  data  adds  to  the  difficulties  involved 
in  predicting  dynamic  response. 
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To  study  the  combined  effects  of  the  new  unsteady  data  and  variable 
inflow,  two  of  the  NH-3A  cases  were  rerun  with  nondistorted  wake  inflow 
distributions.  As  was  found  with  the  original  unsteady  data,  greater 
torsional  moment  oscillations  wore  achieved.  The  calculated  results  were 
very  sensitive  to  small  changes  in  rotor  tip -path -plane  orientation. 

Figure  59  for  case  43  shows  that  the  peak-to-peak  torsional  moment  is  well 
matched  by  the  analysis.  The  advancing  blade  oscillations  are  the  result 
of  calculated  vortex-blade  passage  distances  that  are  much  smaller  than 
those  experienced  in  flight.  The  calculation  was  repeated  with  a  ^-degree 
change  in  tip-path-plane  angle,  and  the  advancing  blade  oscillation  was 
diminished.  The  final  result  is  shown  in  Figure  60.  The  same  sensitivity 
to  rotor  inflow  was  found  for  case  39 •  Inflows  were  calculated  for  the 
measured  tip-path-plane  angle  and  for  a  tip-path -plane  angle  perturbed  by 
1  degree.  The  results  do  not  show  as  good  phase  correlation,  but  the 
sensitivity  of  the  calculation  to  inflow  variations  is  similar  to  that 
shown  for  case  43.  Figure  6l  shows  the  flatwise  and  edgewise  stress 
correlation  obtained  for  case  43  with  the  modified  nondistorted  wake.  In 
this  case,  the  flatwise  vibratory  stress  is  overestimated  but  the  edgewise 
response  is  near  the  level  measured  in  test.  The  nature  of  the  edgewise 
response  lends  support  to  the  feeling  that  proper  modelling  of  wake  effects 
is  important  for  these  flight  conditions. 
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The  determination  of  the  important  effect  of  the  nonuniform  inflow  is 
a  major  result  of  the  NH-3A  correlation  study.  For  these  cases,  small 
tip -path -plane  inclinations  and  fairly  small  rotor  loading  lead  to  small 
wake  skew  angles.  Examination  of  a  group  of  NH-3A  cases  which  had  greater 
negative  tip -path -plane  angles  of  attack  bears  out  the  general  conclusion 
that  high  levels  of  torsional  response  are  related  to  blade -wake  inter¬ 
action.  Induced  velocity  fluctuations  on  the  forward  portion  of  the  rotor 
disk  cause  oscillations  in  blade  pitching  moment  well  before  the  blade 
enters  deep  stall.  Examination  of  the  blade  dynamic  equations  indicates 
that  for  these  cases,  the  aerodynamic  moments  applied  to  the  advancing 
blade  can  result  in  small  torsional  oscillations  over  the  nose  of  the  air¬ 
craft  which  serve  to  define  initial  conditions  for  retreating  blade  res¬ 
ponse.  In  the  third  and  fourth  quadrants  of  rotor  azimuth,  the  blade  acts 
in  torsion  like  a  simple  oscillator  with  very  small  or  even  negative  aero¬ 
dynamic  damping.  Unsteady  aerodynamics  are  needed  to  model  the  growth  rate 
of  the  retreating  blade  oscillation.  The  lack  of  phase  correlation  seen 
in  some  of  tne  torsional  response  calculations  and  the  overprediction  of 
flatwise  and  edgewise  stress  shown  in  Figure  6l  point  out  the  difficulties 
involved  in  comt'ning  a  discrete  wake  analysis  and  a  lifting-line 


aerodynamic  model.  Calculation  of  blade  aerodynamics  at  a  set  of  radial 
locations  can  be  a  poor  representation  of  distributed  aerodynamic  loading 
for  abruptly  varying  inflow.  Representation  of  the  blade  as  a  lifting 
surface  rather  than  a  lifting  line  would  appear  to  be  necessary  in  regions 
where  vortices  are  passing  close  to  blades.  Unsteady  drag  data  are  also 
necessary  for  proper  modelling  of  edgewise  response.  (Steady-state  drag 
data  were  used  in  this  analysis.) 

MODEL  ROTOR  TEST  CORRELATION 

The  analytical  unsteady  model  was  applied  to  the  problem  of  predicting 
torsional  oscillations  obtained  under  controlled  wind  tunnel  conditions. 
Reference  27  reports  the  results  zl  z  aeries  of  tests  performed  with  a 
9 -foot -diameter  articulated  model  rotor.  The  four  untwisted  blades  had 
0.353 -foot  chords,  giving  the  rotor  a  solidity  of  0.1.  The  airfoil  section 
was  NACA  0012.  The  blades  were  dynamically  scaled  replicas  of  full-scale 
blades  having  comparable  mass  and  stiffness  distributions.  At  each  of  a 
series  of  advance  ratios,  blade  collective  pitch  was  incremented  until  high 
retreating  blade  torsional  stress  levels  were  encountered.  For  this  rotor 
system,  the  large  torsional  response  generally  initiated  in  the  third 
quadrant,  increased  in  amplitude  through  the  fourth  quadrant,  and  was 
damped  out  before  azimuth  angle  0=90  degrees  was  reached.  Two  of  the 
more  extreme  conditions  were  chosen  for  analytical  study:  case  51  “H> 
which  had  an  advance  ratio  of  /x=  0.504  and  Ct /ct  =  0.08l,  and  case  68-6, 
which  had  p =  0.294  and  CT/cr  =  0.105 . 

In  the  Normal  Modes  Program,  four  flatwise,  two  edgewise,  and  two 
torsional  elastic  modes  were  used  to  describe  the  blades.  Examination  of 
low  Reynolds  number  aerodynamic  data  indicated  that  lift  and  moment 
stall  angles  are  reduced  by  about  20  percent  from  the  high  Reynolds 
number  data  characteristic  of  full-scale  rotor  operation.  Accordingly, 
the  unsteady  data  were  scaled  to  lift  and  moment  curves  with  stall  angles 
reduced  20  percent.  The  program  was  run  using  measured  cyclic  pitch  and 
shaft  angle.  Attempts  to  increase  collective  pitch  to  attain  measured 
lift  values  were  generally  not  successful.  The  two  conditions  studied 
were  heavily  stalled,  so  increases  in  collective  pitch  did  not  significant¬ 
ly  increase  calculated  lift.  Analytic  lift  was  up  to  15  percent  below  test 
lift  for  nhe  same  collective  pitch  settings. 

When  the  new  unsteady  data  tables  were  used  in  the  prediction  of  model 
rotor  stall  flutter,  too  little  torsional  excitation  was  calculated. 

Figure  62  shows  only  a  small  fourth -quadrant  oscillation  for  case  68-6. 

To  study  the  reason  for  the  small  rate  of  growth  of  the  oscillation,  the 
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two  cases  were  rerun  with  the  original  modified  unsteady  tables.  Figures 
63  and  64  shew  that  the  original  unsteady  aerodynamic  tables  predict  less 
fourth -quadrant  aerodynamic  damping  and  more  nearly  match  the  magnitude 
of  the  stall  flutter  oscillation. 

DISCUSSION  OF  THE  MODEL  ROTOR  CORRELATION 

Unlike  the  NH-3A,  the  model  rotor  seems  to  exhibit  pure  stall  flutter. 
Significant  torsional  excitation  grows  quickly  in  the  fourth  quadrant, 
where  the  blade  angles  of  attack  penetrate  into  stall.  Examination  of  the 
applied  aerodynamic  pitching  moment  and  the  torsional  response  shows 
definite  regions  of  negative  damping  in  the  fourth  quadrant.  A  better 
prediction  of  response  is  obtained  with  the  original  unsteady  data  than 
with  the  new  data  because  the  original  data  predicts  larger  negative 
damping  regions.  This  may  be  due  to  the  low  aspect  ratio  of  the  model 
described  earlier  in  this  report. 

CH-53A  FLIGHT  TEST  CORRELATION 

In  heavily  loaded  flight  conditions,  the  CH-53A  experiences  large 
torsional  oscillations  at  the  blade  torsional  natural  frequency  (approxi¬ 
mately  6. 5 /rev).  These  oscillations  generally  start  over  the  nose  of  the 
aircraft  and  grow  in  amplitude  in  the  third  and  fourth  quadrants  of  rotor 
azimuth.  Case  55  of  the  CH-53A  rotor  loads  report  (Reference  26)  was 
selected  for  analysis.  This  is  a  137”knot  flight  condition  with  C<r/cr  = 
O.O69.  Three  flatwise,  two  edgewise,  and  one  torsional  mode  were  found  to 
be  sufficient  in  describing  the  blade  aeroelastic  response.  This  case  was 
studied  using  test  values  of  shaft  angle  and  cyclic  pitch  while  iterating 
to  the  proper  lift  with  collective  pitch  angle.  Based  on  the  model  rotor 
results,  it  was  decided  to  use  only  the  modified  original  unsteady  data 
for  the  CH-53A  correlation. 

First,  a  series  of  cases  was  run  with  constant  inflow.  As  in  the  case 
of  the  NH-3A,  predicted  torsional  response  was  much  too  small.  Two  more 
constant  inflow  cases  were  run  in  which  the  blade  stall  angle  was  reduced 
by  10  and  20  percent.  Figure  65  shows  the  effect  of  reducing  the  stall 
angle  20  percent.  Deeper  stall  and  higher  levels  of  torsional  excitation 
are  calculated,  but  the  response  does  not  initiate  over  the  nose  of  the 
helicopter  and  the  amplitude  of  oscillation  is  not  matched. 

The  importance  of  nonuniform  rotor  inflow  for  this  case  was  studied 
by  repeating  the  analysis  using  the  variable  inflow  distribution  consistent 
with  a  nondistorted  rotor  wake  and  the  original  blade  stall  angle  values 
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(Table  VIII ).  One  of  the  difficulties  encountered  yith;  the  lifting-line 
analysis  is  that  for  very  close  vqrtex-blad£*passages,  unrealistically 
large  induced  velocities  can  be  calculated.  . Accordingly,  th«  rotor  inflow 
distribution  calculated  from  the  analysis  was  modified  to  /smooth  out  an 
induced  velocity  spike  calculated  at  ip  =  lUO  degrees  and  F  =  0.75*  Figure 
66  shows  the  torsional  signature  resulting  from  this  calculation.  Although 
the  peak-to-peak  torsional  moment  is  not  quite  matched ,  the  phasing  and 
frequency  content  of  the  solution  are  good.  Figure  67  compares  the  root 
torsional  moment  calculated  with  steady  and  unsteady  aerodynamics  for  the 
same  modified  nondistorted  wake  inflow.  With  steady-stdte  aerodynamic 
data,  the  oscillation  begins  similarly  in  the  region  of  tf/  =  180  degrees 
but  damps' out  over  the  retreating  side  of  the  rotor  disk. 

} 

DISCUSSION  OF  THE  CH-53A  CORRELATION 

*  J  ! 

The  analytical  results  for  the  CH-53A  are  similar  to  those  for  the 
NH-3A.  Even  at  this  forward  speed,  a  rotor  tiprpath-plane  angle  of  attack 
of  only  -2  to  -3  degrees  existed.  The  wake  deflection  angle  (angle  between 
the  rotor  disk  and  wake  streamlines)  was  calculated  in  the  Prescribed  Wake 
Analysis  to  be  5.0  degrees.  As  seen  schematically  in  Figure  68,  if  the 
rotor  skew  angle  is  less  than  or  equal  to  the  coning  angle,  tip  vortices 
trailed  from  the  forward  portion  of  the  rotor  difek  can: intersect  inboard 
blade  segments.  Also  shown  in  this  figure  is  the  pattern  of  trailing 
vortices  at  one  instant  of  time  for  case  55.  Several  potentially  strong 
vortex-blade  interactions  over  the  front  half  of  the  rotor  disk  are  shown. 
As  in  the  case  of  the  NH-3A ,  the  observed 1  torsional  response  is  a  hybrid 
made  up  of  wake-induced  effects  :and  retreating  blade  stall  flutter.  Blade 
torsion  is  excited  ovpr  the  nose  of  the  aircraft,  and  the  retreating  blade 
oscillates  in  undamped  motion  at  the  blade  natural  frequency.  The  failure 
to  match  the  amplitude  of  the  test  moment  nay  be  due  to  difficulties  in 
scaling  the  unsteady  pitching  moment  data  to  high  Mach  numbers.  Compared 
to  the  NH-3A  cases,  this  condition  has  lower  forward  speed  and  higher  tip 
speed.  Doth  of  which  increase  retreating  blade  Mach  numbers.  Mach  number 
is  known  to  have  significant  effects  upon  unsteady  characteristics,  and 
correlation  of  high- Mach- number  hysteresis  loops  using  the  present  scaling 

technique  has  been  only  fair. . 

,  •  | 

GENERAL  DISCUSSION  .OF  CORRELATION  . 

The  complexity  of  the  rotor  unsteady  ^hree- dimensional  aerodynamic  . 
environment  suggests  that  complete  understanding  of  blade  torsional 
response  and  control  load  buildup  will  require  further  fundamental  study. 
Blade- wake  interaction,  blade  dynamic  [response  and  dynamic  stall  effects 
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all  contribute  to  the  retreating  blade  torsional  loads.  The  rotor  trail¬ 
ing  wake  was  found  to  be  an  important  source  of  torsional  excitation  not 
only  for  the  NH-3A  compound  but  also  for  the  CH-53A  in  the  pure  helicopter 
mode.  Any  lift  or  propulsive  force  condition  which  has  a  rotor  coning 
angle  greater  than  or  approximately  equal  to  the  wake  skew  angle  can 
result  in  inboard  blade- vortex  intersections.  The  CH-53A  results  indicate 
that  for  the  large  coning  angles  and  small  tip  path  plane  angles  character¬ 
istic  of  high  lift  flight,  blade- vortex  interference  effects  can  become 
significant.  The  effect  of  blade- wake  intersection  on  blade  torsional 
response  is  an  excitation  of  the  blade  first  torsional  mode  by  the  large 
aerodynamic  pitching  moments  resulting  from  vortex- induced  local  stall 
as  the  blade  passes  over  the  nose  of  the  aircraft.  In  response  to  this 
applied  moment  the  retreating  blade  acts  like  a  simple  spring- inertia 
torsional  system  with  small  damping.  The  negative  pitch  damping  associated 
with  dynamic  stall  serves  to  further  increase  the  amplitude  of  retreating 
blade  oscillations. 

Blade  aerodynamics  were  best  described  with  the  original  unsteady 
data  tables  apparently  because  of  the  lower  effective  aspect  ratio  of  the 
more  recently  tested  model  airfoil.  Although  the  analysis  generally 
underpredicted  peak-to-peak  control  loads,  good  correlation  of  the  azimuth¬ 
al  load  distribution  indicates  that  blade  aerodynamic/dynamic  effects  are 
being  reasonably  well  modelled.  For  this  reason  it  is  felt  that  qualita¬ 
tive  predictions  of  control  load  trends  can  be  obtained  with  the  present 
analytical  methods.  However,  continued  correlation  studies  are,  of  course, 
desirable  to  examine  other  parameters  and  to  increase  the  confidence  level 
in  the  analysis.  In  addition,  in  applying  the  analysis,  its  assumptions 
should  always  be  kept  in  mind.  In  particular,  refinements  to  the  blade 
lifting-line  model  used  in  the  analysis  should  be  considered  to  permit 
more  accurate  predictions  of  the  loads  induced  by  a  vertex  passing  close 
to  the  blade. 


CONCLUSIONS 


The  conclusions  listed  below  were  reached  in  the  course  of  this 
investigation.  Although  conclusions  1,  2,  and  3  derive  from  the  current 
data,  it  is  believed  that  these  conclusions  apply  to  the  use  of  general¬ 
ized  data  obtained  from  other  sinusoidal  experiments  as  well.  Con¬ 
clusions  4  and  5  were  reached  upon  performing  correlation  studies  with 
various  rotor  configurations  involving  the  use  of  both  the  new  generalized 
unsteady  data  and  the  original  generalized  unsteady  data  from  Reference  8. 

1.  The  generalized  normal  force  and  moment  data  can  be  used  to 
reconstruct  accurately  the  aerodynamic  response  of  an  airfoil  oscillating 
either  sinusoidally  or  nonsinusoidally. 

2.  The  generalized  moment  data  can  be  used  to  predict  the  stall 
flutter  response  of  a  flexured  airfoil  periodically  driven  into  the  stall 
region  at  low  frequency. 

3.  Pressure  waves  which  may  be  associated  with  dynamic  stall  vortex 
shedding  appear  to  move  rearward  along  the  blade  chord.  Most  of  the 
dynamic  stall  activity  appears  to  be  confined  to  the  airfoil  suction  sur¬ 
face. 

4.  Greatly  improved  correlation  between  measured  rotor  control  loads 
and  those  predicted  by  the  Normal  Modes  Blade  Aeroelastic  Analysis  was 
achieved  relative  to  the  results  obtained  in  Reference  8. 

5.  High- amplitude  oscillatory  control  loads  on  rotor  blades  generally 
result  from  a  combination  of  moments  produced  by  local  stall  induced  by 
blade- vortex  interaction,  stall  flutter,  and  blade  torsional  dynamic 
response.  The  relative  importance  of  these  moment  sources  depends  on  the 
rotor  operating  condition.  Best  results  were  obtained  with  a  smoothed 
version  of  the  original  unsteady  data  tables  and  variable  rotor  inflow. 
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RECOMMENDATIONS 


The  recommendations  listed  below  are  included  here  to  provide  a 
framework  for  the  formulation  of  future  programs  in  the  study  of  dynamic 
stall  aerodynamics  and  its  relation  to  rotor  blade  control  instabilities. 

1.  Analyses  should  be  conducted  to  systematically  define  the  rotor 
design  factors  influencing  the  buildup  in  control  loads. 

2.  Additional  stall  flutter  tests  should  be  conducted  on  a  flexured 
model  with  low  mechanical  damping  driven  at  low  frequency  to  simulate 
cyclic  pitch  changes.  Use  should  be  made  of  the  techniques  developed 
herein  and  a  companion  analysis  performed  to  predict  the  motion. 

3.  Additional  detailed  pressure  measurements  should  be  made  on 
instrumented  models  at  higher  Mach  numbers  to  determine  the  effects  of 
compressibility  on  dynamic  stall. 

4.  The  possibility  of  developing  alternate  methods  for  representing 
rotor  blade  unsteady  aerodynamics  should  be  investigated.  These  methods 
should  strive  to  retain  the  general  accuracy  demonstrated  by  the  a,  A,  B 
method  while  at  the  same  time  avoiding  its  associated  large  data  acquisi¬ 
tion  and  reduction  requirements. 

5.  Unsteady  drag  characteristics  should  be  measured  and  incorporated 
into  the  blade  response  calculations. 

6.  Improved  procedures  for  accurately  predicting  lifts  and  moments 
on  a  blade  as  it  passes  near  a  trailing  vortex  should  be  developed. 

7.  Parametric  stall  flutter  analyses  of  a  two-dimensional  airfoil 
should  be  performed  using  the  analytical  techniques  developed  herein  to 
determine  the  governing  parameters  and  their  effects  on  stall  flutter. 

8.  Pressure  data  obtained  in  this  investigation  should  be  further 
analyzed  to  study  in  detail  the  apparent  pressure  wave  associated  with 
dynamic  stall. 

9.  Additional  stall  flutter  tests  should  be  conducted  with  the 
object  of  visualizing  the  flow  over  the  airfoil  during  the  motion.  This 
could  range  in  sophistication  from  tuft  and  oil  spray  studies  to  high-speed 
Schlieren  motion  pictures. 

38 


10.  Additional  studies  should  be  undertaken  to  further  evaluate  the 
ramp  data  obtained  during  this  program.  Particular  attention  should  be 
paid  to  comparisons  of  force  and  moment  response  over  regions  of  constant 
A. 
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APPENDIX  I 


An  excellent  discussion  of  the  use  of  the  Gaussian  quadrature 
technique  to  integrate  experimentally  measured  chordwise  pressures  is 
given  by  Davis  in  Reference  13.  An  attempt  will  be  made  here  to  summarize 
it  briefly  before  describing  the  application  of  this  technique  to  the 
present  study. 

First,  the  pressure  orifices  must  be  located  at  a  series  of  coordinate 
points  corresponding  to  the  roots  of  a  Legendre  polynomial,  suitably 
normalized  over  the  region  of  integration.  According  to  Davis,  the  method 
has  the  advantage  that  for  n  pressure  stations,  an  exact  integration  is 
possible  if  the  pressure  distribution  can  be  represented  by  a  polynomial 
of  (2n-l)th  degree.  If  the  measured  pressure  distribution  can  only  be 
approximated  by  a  polynomial,  Gauss's  method  will  always  lead  to  a  higher 
order  of  accuracy  for  the  same  number  of  pressure  stations  than  other 
numerical  integration  procedures,  such  as  trapezoidal  or  Simpson's  rule 
integration.  It  is  shown  in  Figure  10  of  Reference  13  that  the  probable 
error  in  computing  the  chordwise  normal  force  using  10  or  more  Gaussian 
stations  is  1  percent  or  less. 

It  was  necessary  to  adapt  the  Gaussian  method  to  apply  it  to  the 
present  experiment.  In  ordinary  usage,  the  formula  for  the  area  under  a 
curve  f(f)  over  the  normalized  interval  0  <  £  <1  is  computed  from  the 
formula 


A--  /  fl^)d$  *£w.fl(.)  (  0<  C<  I )  (10) 

•/  i  =  l  1  ' 

0 

where  the  in  the  right-hand  expression  are  the  Gaussian  weighting  coef¬ 
ficients.  The  N  values  of  ^  sure  located  at  the  zeroes  of  the  Nth  degree 
Legendre  polynomial.  For  10  chordwise  stations,  the  values  of  ^  and 
are  given  in  Table  X.  It  is  seen  from  the  second  column  of  this  table 
that  the  last  f  value  is  located  1.3  percent  from  the  end  of  the  interval, 
and  the  next  to  last  £  value  is  located  6.7  percent  from  the  end  of  the 
interval.  However,  there  was  insufficient  room  in  the  trailing-edge 
region  of  ths  model  rib  to  place  transduce;  <3  aft  of  the  80  percent  chord 
station.  As  a  compromise  solution,  it  was  decided  to  place  the  trans¬ 
ducers  where  there  was  room  ahead  of  the  6u  percent  chord  station,  and  to 
connect  the  three  rearmost  transducers  to  their  measuring  stations  by  means 


of  surface  tunnels  (described  earlier),  but  to  restrict  the  rearmost 
measuring  station  to  be  no  farther  aft  than  the  90  percent  chord  station. 
Thus  it  was  necessary  to  extrapolate  over  the  last  10  percent  of  the  chord 
to  complete  the  chordwise  integrations. 


The  actual  procedure  is  indicated  schematically  in  Figure  69  and  is 
described  below.  Let  rj  be  a  dimensionless  chordwise  coordinate,  0  <  1?  <  1, 
and  let  z<tj  be  some  point  near  the  trailing  edge.  If  Cp(rj)  is  the 
chordwise  pressure  coefficient,  then  the  normal  force  will  be  given  by 
the  integral 


C  (17)  di?  s  /  CArf)  dv  +  /  0.(17) di7 


where  the  right-hand  formulation  represents  the  Gaussian  integration  over 
0<i?<z  and  the  extrapolated  integration  over  z  <7)  <1.  If  we  let 


rr-  z( 


Equation  (ll)  becomes 


+/ c, 


,(17)  ^17 


Replace  the  first  integral  of  Equation  (13)  by  its  Gaussian  equivalent 
from  Equation  (10), 


Cn=zEwiCp(Ci)  +1  CAr ?)dV 


The  rearmost  orifice  was  chosen  to  be  located  at  the  90  percent  chore, 
station  (see  Table  X,  right-hand  column),  so  =  0. 98695  corresponds  to 
17=  O.9O.  Hence,  from  Equation  (12),  we  can  solve  for  z  to  obtain  z  = 
0.912.  Equation  (l4)  then  becomes 


cn*  0.912 


(15) 


t 


+/c 

.912 


To  evaluate  the  last  integral,  assume  a  parabolic  form  to  Cp(i7), 


Cp(i?) 


(16) 


and  require  that  the  curve  pass  through  the  coordinates  (179,  Cp^),  (t^q, 
Cp10)  and  (1,  0).  The  final  result  for  this  integral  will  then  be 


/ 


Cp(i))dT)  *  CA(I-  zl  +  i-  Cb(I-Z2>  +  Cj.U-25) 


z 


(17) 


where 


A  = 


179)(l-T?tO)(l_179) 


(18) 
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and  where 


Z  1  0.9119 

*  «  0.850374 
'9 

V  0>9° 

The  moment  coefficient  about  the  25  percent  chord  station  is  obtained 
from  these  equations  by  replacing  Cp(^)  and  Cp(7j)  by 

F(Ci>*(0.25-£j)Cp(£i>  (20) 

and 

nv)'  (0.25-t,)Cp(^)  (21) 

respectively. 


47 


APPENDIX  II 


OPTIMUM  RECORDING  CONDITIONS  FOR  WISARD 


Fourier  spectral  methods  were  not  directly  used  in  the  analysis  of 
the  unsteady  data.  However,  the  ability  of  the  data  system  to  extract 
information  with  sufficient  frequency  resolution  to  ensure  accuracy  and 
repeatability  was  dependent  on  many  of  the  same  parameters  that  govern  the 
resolving  capabilities  of  the  spectral  analysis  technique  (Reference  30). 
Hence,  the  parameter  values  selected  for  the  data  recording  mode  were  the 
same  as  those  that  would  be  required  for  a  spectral  analysis. 

In  this  brief  analysis,  we  seek  to  establish  the  criteria  which  were 
used  to  select  the  tape  speed  at  which  data  were  recorded  (SR  in  in. /sec) 
and  total  record  time  (T  in  sec).  To  do  this,  we  must  first  note  that  the 
absolute  sampling  rate  of  WISARD  during  data  playback  (which  is  the  rate 
at  which  data  were  digitized  and  stored  on  digital  tape)  was  15,000 
samples/sec.  This  rate  was  divided  by  the  number  of  channels  being 
sampled,  Nq,  to  obtain  the  absolute  sampling  rate  pier  channel  during  play¬ 
back, 


AT, 


ABS 


15,000 

N, 


(  samples  /  sec ) 


(22) 


where  AT^gq  is  the  absolute  time  between  samples. 


In  general  during  this  test  program,  the  recording  speed,  SR,  was 
greater  than  the  playback  speed,  Sp.  In  this  case,  the  effective  sampling 
rate  of  the  original  data,  within  the  time  frame  of  the  original  data, 
was  greater  by  the  ratio  SR/Sp,  or 


h 


R 


_ ! _  x  Is 

atabs  sp 


(23) 


where  hR  is  the  effective  sampling  rate  and  AT  is  the  time  between  samples 
in  the  original  time  domain.  Hence,  from  Equations  (22)  and  (23), 
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ggrar  Ena— 


.  .  I  .15.000 

*  AT,  Nr 


(  samples  /  sec) 


In  many  cases,  the  core  storage  capacity  of  the  computing  machine 
was  a  limiting  factof.  In  the  spectral  analysis  procedure  outlined  in 
Reference  30,  a  specific, parameter,  m,  is  denoted  as  the  total  number  of 
samples  to  be  processed  per  channel.  This  is  most  conveniently  chosen  to 
be  some  integral  power  of  2,  the  most  useful  of  which  are 


512  =  2' 


samples  /channel 


1024-  2  u 
2048=  2  11 
•-J2 

4096=  2 

„I3 

8192  2 


Thus  the  minimum  recording  time  to  obtain  m  samples  per  channel  is 
obtained  by  dividing  Equation  (25)  by  Equation  (24),  or 


T=  _m_  = _ V- -c  ...  ... 

hR  15,000  (SR/Sp) 


( sec ) 


In  Reference  30,  it  is  pointed  Out  that  an  important  parameter  in  the 
analysis  is  the  minimum  number  of  samples  per  cycle  of  the  highest  har¬ 
monic  of  interest;  this  parameter  is  denoted  as  L  in  the  present  work.  A 
criterion  is  established  in  Reference  1 30  that  when  L  =  2  samples/cycle,  a 
maximum  frequency  can  be  determined  from  Equation  (24), 


hR  -  hR 


(cps) 


which  sets  an  upper  lipit  for  meaningful  results.  This  is  the  highest 
frequency  attainable  before  aliasing  takes  pla^e,  in  which  higher  frequency 
information  (at  f  >fmax)  is  folded  into  the  low-frequency  domain 


P 


(see  Reference  30  for  more  details  on  aliasing).  Use  of  a  minimum  value 
of  L  =  4  samples/cycle  in  the  present  analysis  and  a  low -pass  filter  to 
block  the  high-frequency  noise  ensured  the  accuracy  of  the  data  to  well 
beyond  the  highest  harmonic  specified  as  fmaY . 

A  final  parameter  of  interest  is  the  bandwidth  which  defines  the 
frequency  resolution  of  the  computation.  This  is  defined  in  Reference  30 
as  twice  the  maximum  frequency  of  interest  divided  by  the  number  of 
samples  being  analyzed,  or 


B 


2  f, 


max 


w 


m 


(28) 


When  Equations  (24)  and  (27)  (with  L  =  4)  are  substituted  into  Equation 
(28),  the  result  is 


B  :  7500  /  Sr  ) 


'w 


mN, 


(29) 


The  number  of  channels  in  this  application  was  Nc  =  12,  and  the 
number  of  samples  per  channel  chosen  for  analysis  was  m  =  4096.  With  the 
introduction  of  these  numbers,.  Equations  (26),  (28),  and  (29)  can  be 
rearranged  to  yield 


!a_  = 


f 


max 


312.5 


(30) 


Bw=  0.15259 


(3D 


T  = 


3.2768 

(vv 


(32) 
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These  equations  were  used  to  construct  Table  III.  The  highest  harmonic  of 
interest,  fmav .  was  chosen  to  be  10  times  the  test  frequency,  and  Sp/Sp 
was  computed  from  Equation  (30).  Use  of  a  fixed  value,  Sp  =  1.875  in. /sec, 
fixed  the  calculated  value  of  Sp.  However,  the  tape  recorder  could  operate 
only  at  fixed  multiples  of  the  playback  speed,  so  the  next  highest  value 
of  Sp  was  chosen  and  the  actual  speed  ratio  Sp/Sp  was  computed.  This  was 
then  introduced  into  Equations  (31)  and  (32)  to  compute  the  frequency 
resolution  and  record  time.  Finally,  the  actual  number  of  samples  per 
cycle  in  the  10th  harmonic,  was  obtained  by  solving  Equation  (27) 

for  L,  substituting  Equation  (24)  for  hR,  and  using  Nc  =  12  and  the  actual 
tape  speed  ratio, 


-  1250 


(t) 

\  'ACT 
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The  single -degree -of -freedom  system  consisting  of  the  NACA  0012  air¬ 
foil  mounted  in  the  wind  tunnel  on  a  torsional  spring  was  analyzed  using 
the  a,  A,  B  tabulation  of  unsteady  moment.  In  the  test,  the  airfoil  was 
held  in  the  tunnel  sidewall  plates  by  thrust  bearings  which  introduced  a 
significant  amount  of  damping.  As  mentioned  in  the  text,  the  airfoil  was 
excited  sinusoidally. 

The  equation  of  motion  of  the  mechanical  system  described  is  the  same 
as  that  for  a  single -degree-of -freedom  torsional  forced  vibration  problem 
with  damping  but  also  with  an  added  external  moment  resulting  from  the 
airflow  over  the  airfoil.  The  equation  of  motion  may  be  written  as 


la  +  Ca  +  K(a-aM>  =  M(t, a, a, «!,...)+  T0cos  (34) 


where 


I  =  airfoil  mass  moment  of  inertia  in  pitch 
C  =  equivalent  viscous  damping 
K  =  torsional  spring  constant 
M(  )  =  aerodynamic  moment 

T0  =  magnitude  of  the  applied  torque 
<u  =  frequency  of  applied  moment 
t  =  time 

Of  =  instantaneous  angle  of  attack 

=  mean  angle  of  attack  of  the  oscillation 
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The  primary  purpose  of  this  program  was  to  test  the  ability  of  the 
sinusoidally  generated  unsteady  force  and  moment  data  to  predict  non- 
sinusoidal  unsteady  force  and  moment.  Accordingly,  the  aerodynamic  moment 
used  in  the  differential  equation  (Equation  (3*+))  was  assumed  to  be  a 
function  of  angle  of  attack  and  its  first  two  time  derivatives  only,  and 
it  was  expressed  as  M(a,  A,  B).  Equation  (3*0  may  be  divided  through  by 
2pV2b2  to  convert  the  moment  to  the  moment  coefficient,  with  the  result 


2pV2  b2 


a  + 


f  la-v 


=  CmU,A,B>  + 


T0  COS  cut 
ZpV2  b* 


(35) 


If  we  define  the  quantities 


I 


a 


I 

2pb4 


(36) 


C_  .. 

I  IT 


(37) 


K_ 

I 


(38) 


where  o>ais  the  torsional  natural  frequen  y  and  8  is  the  logarithmic 
decrement  of  the  damped  oscillation,  then  Equation  (35)  may  be  rewritten 
as 


I 


a 


w„  (  a-a  ) 
a  m 


CM(a,AtB) 


T0  cos  (vt 

+  2  pVZ  b2 


(39) 


Dividing  through  by  Ia(b/v)^  yields 


(a-aM)  = 


Cm(a,A,B)  T0  cos  qjt 
Ia(b/V)2  +  I 


(40) 
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This  second-order  differential  equation  can  be  converted  to  two  first-order 
equations  for  solution  on  the  computer  by  assuming  that 


a  ■-  $ 


(41) 


Now  Equation  (4o)  may  be  written  as 


a  -a 


Cm(a,A,B) 
Ia(  b/y)2 


+ 


T0  COS  cut 
I 


(42) 


Equations  (4l)  and  (42)  are  a  pair  of  first-order  differential  equations 
suitable  for  solution  on  the  computer.  It  should  be  noted  that  since 
A  =  ba/v  and  B  =  b2/j/v2,  the  moment  coefficient  is  a  function  of  the  para¬ 
meters  being  solved  for.  The  assumption  has  been  made  in  the  solution 
that  if  a  sufficiently  short  time  step  is  taken,  the  moment  coefficient  is 
adequately  represented  by  using  past  values  of  A  and  B  when  interpolating 
Table  V  for  C^. 

The  initial  conditions  are  taken  at  the  peak  of  a  cycle  of  motion  as 


an  =  a(0)  =  a,  +  a 

U  M 


(43) 


/305  /3(0)=  0 


(44) 


The  solution  was  started  with  the  best  estimate  of  the  actual  Cm,  namely, 


m 


C.  ( 


m 


a,  + 


0,  - 


k2  2  _ 
d  ui  a 


/V2) 


(45) 


•4 


By  substituting  Equations  (48)  and  (50)  into  Equations  (46)  and  (47) 
together  with  the  initial  conditions,  the  quantities  aQ,  (3Q,  a and 
are  completely  defined. 


The  numerical  integration  routine  can  now  be  started  to  obtain  the 
complete  solution.  A  simple  trapezoidal  formula  (see  Reference  31)  was 
used  with  the  following  recursion  relationships : 

a  =  <*„  +  +  £„)  (n) 
n+i  n  2  ^n+i  V5J-; 


^n  +  ,=  ^n  + 


At 

2 


1 


(52) 


where  was  obtained  from  the  original  differential  equation  (Equation 
(42))  with  continually  updated  values  of  Cm. 

The  numerical  integration  was  continued  until  two  consecutive  cycles 
of  motion  had  converged  to  within  an  angle  of  attack  discrepancy  of  0.5 
percent.  Original  computations  were  made  to  a  tolerance  of  0.1  percent, 
but  the  added  accuracy  occasionally  caused  convergence  difficulties.  These 
problems  were  most  probably  caused  because  the  tabulated  moment 
coefficients  had  small  fluctuations  and  were  not  completely  smooth  mathe¬ 
matical  functions .  The  larger  tolerance  amounted  to  a  maximum  of  only 
0.1  degree,  which  was  roughly  the  plotting  accuracy  of  the  data  and 
therefore  was  not  significant. 
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TABLE  I.  TRANSDUCER  SENSITIVITIES 


Position 

Number 

Percent  Chord 

Behind  L.E. 

Sensitivity 

(mv/psi) 

1 

1.19 

1.81 

2 

6.15 

2.06 

3 

14.62 

1.23 

4 

25.83 

1.385 

5 

38.81 

2.385 

6 

52.38 

I.85 

7 

65.36 

1.58 

8 

76.57 

1.295 

9 

85.03 

O.98 

10 


90.00 


1.76 


Iter. 

No. 

Can  Qescrli  tlon 
and  Angular 
Amplitude 

Type  of 
Coupler 

Primary* 

drive 

Frequency 

(cps) 

Mean  Angle* 
of  Attack 

(deg) 

1 

none 

rigid 

0  (steady  state) 

0,3,7,9,10,11, 
12, 13, 14, 16, 

lb, 20,22 

1  - 

sinusoidal,  *B 

rigid 

12.5,  31,  50,  75 
100 

3,0,9,11,12, 

It*, 16, lb 

i 

sinusoidal,  1 1 

rigid 

120 

3 

4 

sinusoidal,  *S 

rigid 

50,75 

11  (single 
surface  mean. ) 

e 

sinust  iial,  t  3 

rigid 

10,25,50,75 

14 

i 

tielice;  tar,  *  b 

rigid 

10,16. 7,33 

6,11,14 

7 

forward  ran;  , 

*fc 

rigi  1 

7.5,10,15,20 

o,ll, la 

fc 

backvari  ran;  , 

±  * 

rigid 

7.5,10,15,20 

*  ,11,14 

u 

sinusoidal,  ±t 

flexible 
( 3  s'  rings' 

7.5,10,15 

,11,1c 

lo 

no  ne 

fler.it  le 
t  3  s t  rings 

0  < free  vibrati  n 

0,11, It  *  <» 
for  ."ax. 
-scillati.n 

11 

:'.e 

flexible 
(3  springs 

intuitive  l..a  tin,;, 

An  3 

6,11,16 

•  ”  ‘.e: 

3acr.  iter,  rer resents  test 
angle  sf  attacr.. 

points  taken  at  all 

possible  combinations  f 

fre  pjency  an  1 
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III.  TAPE  RECORDING  PARAMETERS  FOR  WISARD 
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TABLE  VII.  PARAMETER  VALUES  FOR  NUMERICAL  SOLUTION 


I. 

1 


I 


I 

= 

0.000529  (ft«d»b-sec2/rad)/ft 

b 

= 

0.2x3  ft 

V 

= 

390  ft/sec 

P 

= 

0.0022  lb  sec2 /ft1* 

To 

= 

=  0.739  un2  x  10"' ^  ft-lb/ft 

6 

= 

various  values  from  0.25  to  1.00 

u 

= 

47.1,  62.8,  and  92.4  rad/sec  \ 

Un 

= 

260.5,  326,  and  471  rad/sec  v 

au 

= 

6,  31,  and  16  deg  1 

a 

= 

8  deg 

At 

= 

0.002 u  =  3.19w  x  10“^  sec 

2  TT 

Tolerance 

= 

|“in+l  '  a  ini  <  °-°°5 

ain+l 


All  Combinations 


TABLE  VIII.  GCAL  DIG 

PARAMETERS  USED  FOR  MOEEL  AIRFOIL  DATA 

Mach 

Number 

M 

Normal  Force 
Stall  Angle, 

a  sn  (deg) 

Normal  Force 

Coefficient 
at  Stall, 

*n 

Pitching  Moment 
Stall  Angle, 

ofsm  (deg) 

0 

11.0 

1.17 

13.5 

0.2 

11.0 

1.17 

14.0 

0.3 

11.0 

1.248 

14.0 

0.4 

9.0 

1.078 

12.4 

0.5 

6.8 

0.845 

9.0 

0.6 

5.6 

0.740 

8.0 

0.7 

4.0 

0.528 

7.0 

0.8o 

2.0 

0.418 

7.0 

TABLE  IX.  SUMMARY  OF  TEST  CASES  USED  FOR  CORRELATION  STUD* 


TABLE  X.  COORDINATES  AND  WEIGHTING  COEFFICIENTS 

FOR  10  POINT  GAUSSIAN  QUADRATURE 

Station 

i 

Normalized 

Coordinate 

«i 

Weighting 

Coefficient 

Wi 

Actual  Chordwise 
Location, 
Percent  Chord 

1 

0.01305 

0.033336 

1.19 

2 

0.06747 

0.074726 

6.15 

3 

0.16030 

0.109543 

14.62 

4 

0.28330 

0.134634 

25.83  1 

5 

0.42557 

0.147762 

38.81 

6 

0.57443 

0.147762 

52.38 

7 

0.71670 

0.134634 

65.36 

8 

0.83970 

0.109543 

76.57 

9 

0.93? 53 

0.074726 

85.04 

.10 

0.98695 

0.033336 

90.00 

L^Uyu 


ADJUSTABLE 

QUADRANT 


LOCATION  OF 
ISOLATED  AIRFOIL 


PIVOT  AXIS  FOR 
SIDEWALL  QUADRANT 


Figure  2.  UAC  Two-DimenBional  High-Speed  Cascade  Wind  Tunnel 
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INTERIOR  STRUCTURE  OF  AIRFOIL 


Figure  3*  Airfoil  Model  Rib  and  Spar  Assembly 
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MINIATURE  DIFFERENTIAL  PRESSURE  TRANSDUCER 


i 


TRANSVERSE  SECTION  VIEW  THROUGH  LARGE  RIB,  LOOKING  AFT 
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Figure  7.  Miniature  Pressure  Transducer  Installation. 


'igure  o.  Instrumen' 
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BEARING 


Figure  9.  Model  Drive  System  -  Schematic 
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DETAIL  OF  INTERIOR  OF  TRANSMISSION 
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Figure  110.  Capacitive  Angular  Transducer  -  Schematic 


PRESSURE  DIFFERENCE 


O  STATIC 


Time  History  for  a  Single  Cycle  in  Separated  Flow. 


For  each  a  ,  v 


•  Choose  a  value  of  a 


•  Choose  a  value  of  B 


— Select  anaM  such  that  la  —  aMl  <  <2 

•  Each  aM  determines  a  value  of  f  from  Eq.  (38) 
and  A  from  Eq.  (39)  ,  both  from  Ref.  (8) 

•  Enter  appropriate  cross-plot  of  Cmvs  f  ;  read 
Cmat  calculated  value  of  f 

•  Plot  Cm  vs.  calculated  value  of  A 


■ - •  Return  for  new  value  of  aM 

•  Return  for  new  value  of  B 


•  Return  for  new  value  of  a 


Figure  IT-  Flow  Chart  for  Data  Conversion  to  a,  A,  B  Form. 
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Figure  l8.  Steady-State  Normal  Force  vs.  Angle  of  Attack. 


0.05 


NONDIMEN  SIONAL  TIME 


Figure  21.  Variation  of  Angle  of  Attack,  Angular  Velocity  Parameter,  and 

Angular  Acceleration  Parameter  With  Time  for  Typical  Ramp  Motion 


BACKWARD  RUNNING  RAMP  ?  F  -  20  CPS 


Figure  2k. 


ANGLE  OF  ATTACK  INCREASING  FROM 
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Figure  25.  Instantaneous  Chordvise  Pressure  Distribution  on  Pressure 
Surface  During  a  Sinusoidal  Oscillation. 
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ANGULAR  ACCELERATION  PARAMETER,  B  =  b  u/V 


CURVE 


MOTION 


NORMAL  FORCE  COEFFICIENT 


PITCHING  MOMENT 
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A  RECONSTRUCTED  LOOP 
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ANGLE  OF  ATTACK,  «.  DEG 


Figure  29.  Comparison  of  Original  and  Reconstructed  Pitching  Moment 
Hysteresis  Loops  for  Sinusoidal  Motion  at  Various 
Frequencies  and  Mean  Angles  of  Attack  for  a  =  ±8  Degrees. 


PITCHING  MOMENT  COEFFICIENT.  Cm  NORMAL  FORCE  COEFFICIENT 


12.6  CPS 


ore  30-  Comparison  of  Original  and  Reconstructed  Normal  Force  and  Pitching 
Moment  Hysteresis  Loops  for  ±3  Degrees  Sinusoidal  Motion. 
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PITCHING  MOMENT  COEFFICIENT,  C_  NORMAL  FORCE  COEFFICIENT 


COEFFICIENT 


PREDICTED  (BASED  ON  SINUSOIDAL  DATA 
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U  -  47.1  RAD/SEC  (f  -  7.5  CPS) 


v\  NOTE.  06.0n,  AND  014 

\\  REFER  TO  ORIGINS 

'\  FOR  6,  11,  AND  16  DEC 

\\  MEAN  ANCLE  OF  ATTACK 


- ANALYTICAL  RESULTS 
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Figure  kO.  Comparison  of  Experimental  and  Analytical  Time  Histories  for 
Flexible  Airfoil  With  f a  =  52  CPS. 
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DEC 


<*>  -  94.2  RAD/SEC  (f-  15  CPS) 


DEG 


SMOOTHED  DATA 
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NORMAL  FOF  CE  COEFFICIENT,  C„  PITCHING  MOMENT  COEFFICIENT 


NORMAL  FORCE  COEFFICIENT,  C 


0  4  8  12  16  20  24  28 

ANGLE  OF  ATTACK,  «  ,  DEG 
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Figure  b6.  NACA  0012  Normal  Force  Characteristics. 
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NALYS1S  USING  ORIGINAL  MODIFIED  UNSTEADY  AERODYNAMICS 


measured  values 


MEASURED  VALUES:  C  t/ct  =  0.058,  u  =  0.474 


Figure  56.  Torsional  Response  Correlation  Using  Steady  and  Unsteady  Aerodynamics  With  Constant 
Inflow  -  NH-3A  Case  43. 
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Figure  57-  Torsional  Response  Correlation  Using  Steady  and  Unsteady  Aerodynamics  With  Constant 
Inflow  -  NH-3A  Case  72  -  6o  Degree  Angle  of  Bank,  Right  Turn. 
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Figure  58. 


MEASURED  VALUES 


MEASURED  VALUES:  C  T /a  =  0jD58,  ti  =  0.474 


Figure  6l  -  Continued 


ANALYSIS  USING  NEW  UNSTEADY  AERODYNAMIC 


Figure  62.  Torsional  Response  Correlation  Using  New  Unsteady  Aerodynamics  -  Model  Rotor  Case  68-6- 


UES:  6  75  =  11  (l  =  0.294 
EST  CT/ff  =  0.105 

NALYSIS  USING  MODIFIED  ORIGINAL  UNSTEADY 


220  260  300 


MEASURED  VALUES 


AZIMUTH,  0  .  DEC 

'orsional  Response  Correlation  Using  Original  Unsteady  Aerodynamics 


± 


MEASURED  VALUES:  CT/<r=  0.069,  fi  •=  0.324 


Figure  66.  Torsional  Response  Correlation  Using  Original  Unsteady  Aerodynamics  and  Modified 
Nondistorted  Wake  Inflcrw  -  CH-53A  Case  55. 


MEASURED  VALUES:  CT  a  =  0.069,  u  =  0.324 


Figure  67.  Torsional  Response  Calculated  Using  Steady  and  Unsteady  Aerodynamics  and  Modified 
Nondistorted  Wake  Inflow  -  CH-53A  Case  55. 


Wake  Geometry  -  CH-53A  Case  55. 


